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The diffraction of electromagnetic radiation by a hole 
small compared with the wave-length is treated theoreti- 
cally. A complete solution is found satisfying Maxwell's 
equations and the boundary conditions everywhere (Sec- 


tion 4). The solution holds for a circular hole in a perfectly 


conducting plane screen, but it is believed that the method 
will be applicable to much more general problems (Sec- 
tion 8). The method is based on the use of fictitious mag- 
netic charges and currents in the diffracting hole which has 
the advantage of automatically satisfying the boundary 
conditions on the conducting screen. The charges and 
currents are adjusted so as to give the correct tangential 
magnetic, and normal electric, field in the hole. The result 
(Section 5) is completely different from that of Kirchhoff's 


method, giving for the diffracted electric and magnetic 
field values which are smaller in the ratio (radius of the 
hole/wave-length) (Section 6). The diffracted field can be 
considered as caused by a magnetic moment in the plane of 
the hole, and an electric moment perpendicular to it 
(Section 6). The theory is applied to the problem of mutual 
excitation of cavities coupled by small holes (Section 9). 
This leads to equations very similar to those for ordinary 
coupled circuits. The phase and amplitude relations of 
two coupled cavities are not uniquely determined, but 
there are two modes of oscillation, of slightly different 
frequency, for which these relations are opposite (Section 
10). The problem of stepping up the excitation from one 
cavity to another is treated (Section 11). 


1. THE PROBLEM 


N microwave work it is often important to 
know the effect of a small hole in a cavity 
upon the oscillation of that cavity. For instance, 
two cavities may be coupled by a small hole in 
their common boundary (Fig. 1); in this case, 
we wish to know the characteristic frequencies 
and the phase relations for the oscillations of the 
coupled system. Or a hole in a cavity may serve 
the purpose of getting radiation out of it; then 
we want to calculate the amount and the spatial 
distribution of the emitted radiation. Another 
similar problem would be to calculate the effect 
of a small gap in a wave guide upon the propaga- 
tion of waves along that guide. 
A less practical problem but probably the 
simplest one of the same type, is the diffraction 
of electromagnetic waves by a small hole in an 


infinite plane conducting screen. This is the 
problem which we are going to solve first 


Fic. 1. Two cavities, a and 8, coupled by a small hole. 
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(Sections 3-7) ; the application of the solution to 
the practical problems mentioned is then rather 
straightforward (Sections 9-11). 

The available theoretical methods are entirely 
inadequate for the treatment of our problem. 
In the usual Kirchhoff method, the diffracted 


‘ field is expressed in terms of the incident field 


in the hole. However, the Kirchhoff solution does 
not satisfy the boundary conditions, viz., it does 
not give zero tangential component of the 
electric field on the screen. In most textbooks, 
the pious hope is expressed that Kirchhoff's 
method will give at least the first term of a 


‘ convergent series. This is probably true for the 


diffraction by an opening, large compared with 
the wave-length, because then the diffracted 
field will be relatively small on the screen, thus 
“almost” fulfilling the boundary conditions. But 
it is certainly not true for a small hole; in fact, 
our exact solution of the problem will turn out 
to be entirely different from Kirchhoff's. The 
failure of Kirchhoff’s theory will be demonstrated 
mathematically in Section 2. 

Kirchhoff’s method has the additional defect 
of being a scalar theory while the electromagnetic 
field is essentially vectorial. This shortcoming 
has been remedied in the last forty years by a 
number of writers; a very good account of the 
vector equivalent of Kirchhoff’s theory is given 
in Stratton’s book,' to which we shall frequently 


refer in this paper. The vectorial theory ensures 


the fulfillment of the divergence conditions, 
div E=div H=0; i.e., it gives transverse waves 
in the wave zone which would not necessarily 
be the case in Kirchhoff’s scalar formulation. 
However, the vector formulation in no way 
improves the situation regarding the fulfillment 
of the boundary conditions on the conducting 
screen. 

The only rigorous solution of any diffraction 
problem known to me is Sommerfeld’s solution? 
of the diffraction by a conducting semi-infinite 
plane, or by a wedge. As is well known, this 
solution is rather complicated although the 
problem is the simplest imaginable, being two- 
rather than three-dimensional. It appears hope- 


1J. A. Stratton, Electromagnetic Theory 
Book Company, Inc., New York, 1941), p. 464 ff 
2A. Sommerfeld, Riemann-Weber’s Differentialgleich- 


ungen der Physik, seventh edition, p. 433. 
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less to look for a rigorous solution of our problem 
along the lines of Sommerfeld’s solution. A slight 
similarity between his and our method will be 
pointed out later (Section 6d). 

The main simplifying assumption we are going 
to use is that the hole is small compared with 
the wave-length. This means that the incident 
electromagnetic field is almost constant over the 
hole. We believe, however, that our method can 
be generalized to holes of larger size (Section 8). 


2. THE FAILURE OF KIRCHHOFF’S THEORY 


The failure is most easily seen in the simplest 
formulation of the theory, viz., Kirchhoff’s own, 
scalar formulation. Let u be a scalar wave 
function satisfying the wave equation 


Vut+k*u=0. (1) 


We may identify u, e.g., with one of the com- 
ponents of the electric field, say E,. The deriva. 
tives of u will then be connected with the mag- 
netic field. If the conducting screen is placed at 
x=0, the boundary condition is 


u=0 at x=0. (2) 


Let electromagnetic waves ‘come in from the 
left (x <0), the corresponding wave function may 
be uo. Then according to Green’s theorem, the 
wave function at any point r on the right of 
the screen is 


Ouo 


with the Green’s function 
o(r) =e*"/r. (3a) 


The integration variables are y’, 2’. The coordi- 
nates r of the “field point”’ and x’ (=0) are kept 
constant. The solution (3) is rigorous if the 
integral is extended over the entire surface x’ =0, 
with the correct values for 0(0, y’, 2’). However, 
the distribution of uw» on the plane x’ =0 is not 
known. Kirchhoff’s method consists now in 
putting %9=0uo/dx’ =0 everywhere on the screen 
outside the hole, and replacing uo by the field 
of the incident wave in the hole. 

Suppose now the hole is very small compared 
with the wave-length, then u» may be assumed to 
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be constant over the hole and we obtain, re- 

membering that ¢ is a symmetrical function of 
r and r’, 

duo d¢(r) 

dx’ 0 


(4) 


where A is the area of the hole. We can now 
make two alternative assumptions, viz.: 

(a) uo in the hole has the value given by the 
unperturbed incident wave. Then both uo and 
du,/dx’ are different from zero. If we have an 
incident plane wave traveling at an angle 3 with 
the x direction, we have 


(4a) 


(b) We may take into account the reflection of 
the incident wave by the conducting screen and 
put % equal to incident plus reflected wave. 
Then we have in the hole %)=0 while the value 
of 8uo/dx’ is doubled. 

Now consider the expression (4) on the screen 
(x=0). Then it follows immediately from (3a) 
that d¢/8x=0 so that the second term in (4) 
satisfied the boundary condition (2). The first 
term in (4), however, does not vanish by any 
means. Unfortunately, it is just this first term 
which remains present when assumption (0) is 
made. Therefore, for either of the assumptions 
(a) and (6), the boundary condition (2) is 
violated on the “right-hand”’ (back) side of the 
screen. This would only be different if we could 
represent the field by the second term in (4) 
alone. 

The vectorial theory has the same defects. 
If there are no currents or charges in the half- 
space x >0, the electric field at a point 7 in that 
half-space is given by Eq. (19), p.466, in Stratton’s 
book, which is 


=tkuy cos 


E( do(ikn X H(r’ 
n=—f [iknXH(r')o 


—(nXE(r')) Xgrad g—n-E(r’) grad ¢]. (5) 


In this is a unit vector in the direction of the 


inward normal to the surface, i.e., in our case 
the positive x direction. The first integral in 
Stratton’s formula is omitted because there are 
assumed to be no charges or currents. The nota- 
tion is slightly changed (1) because we use 
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Gaussian units so that Stratton’s wy is replaced 
by k, (2) because we use the gradient with respect 
to the coordinates of the field point r (point at 
which the field is to be calculated) rather than 
of the source point r’ (point on the surface where 
the field is given), (3) because we have inter- 
changed the notations r and r’, and (4) because 
we have reversed the sign of n. 

If the integral is extended over the entire 
surface x’=0, Eq. (5) is of course correct. How- 
ever, according to Kirchhoff’s method, the in- 
tegral is only extended over the hole, and for EZ 
and H we insert the field of the wave incident 
from the left, i.e., Eo, Ho. These quantities may 
be again considered‘ as constant over the hole 
which permits integration as in Eq. (4). 

Now let us again consider the electric field at 
a point r on the screen. Then grad ¢ will be a 
vector in the plane of the screen, and the same 
is true of mXH and of XE. Therefore the first 
and last term in (5) give tangential components 
of £, thus violating the boundary condition 
Evan=0. Only the second term satisfies the 
boundary condition, giving an electric field 
normal to the surface. Again, this second term 
is the only one which drops out if we include the 
reflected wave in calculating E(r’) and H(r’) in 
the hole. 


3. MATHEMATICAL FORMULATION 


Let Ho, Eo be the field on the left-hand side 
of the screen if there is no hole. This field fulfills 
the boundary condition for x=0 (plane of the 
screen) 

Eo tan =0 (6) 
which may also be written in the form 
nXE,=0. 
This makes automatically 


Hy-n=0. 


(6a) 


(6b) 


Ho tan, and Eon, are in general different from 
zero. On the right-hand side of the screen, the 
field in zero approximation vanishes identically. 
Then the zero approximation field satisfies the 
boundary conditions everywhere on the screen 
but not in the hole: In the hole, H and the 
normal component of E are discontinuous. 
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We write now the actual field 
for x<0, 
for x>0, 
and similarly for the electric field. Then we have 
the boundary conditions: 
Ei tan=E2 tan in the hole, (8) 


Ei tan=Es2 tan=0 for x=0 outside the hole, (8a) 


A tan — Hy ten = He tan in the hole. (8b) 


The boundary conditions for the normal com- 
ponents are automatically. fulfilled if those for 
the tangential components are satisfied. 

It is easily seen that all boundary conditions 
for and Hy are satisfied if satisfies the 


- boundary condition (8a) and if we put, for any 


x>0, and any y and z: 
y, 2) = Ex, (x, (9) 
Ai,(—x, — Hy, (x, 2), (9a) 


and correspondingly for the z components, Max- 

well’s equations are consistent with (9), (9a) and 
will further make 

2) = H2,(x, 2), (9b) 

y, 2) = —E2z(x, y, 2). (9c) 


Inserting (9a) into (8b) we find that Hz must 
satisfy the boundary condition 


tan=4$Ho tan in the hole. (10) 
Likewise, from (9c) we get the similar condition 
in the hole. (10a) 


Hy and Eo,, may be considered as known. 

The problem is then to calculate the field 
E:, Hz subject to the boundary conditions (8a), 
(10), and (10a). These conditions are valid 
irrespective of size and shape of the hole. How- 
ever, in the following, we shall assume the hole 
to be small compared with the wave-length so 
that Hoy, Hoz, and Eo, may be considered as 
constant over the hole, and we shall take the 
shape of the hole to be circular, of radius a. 


4. SOLUTION 


- As we have seen in Section 2, only one of the 
terms in (5) leads to an acceptable solution, 


vig., the second. Such a term would be produced 
by a distribution of ‘‘magnetic currents” over 
the hole (cf. Stratton’). We shall therefore 
assume a distribution of magnetic currents in the 
plane of the hole, but instead of assuming the 
current density to be proportional to XE, we 
shall determine it so as to satisfy the boundary 
conditions (10), (10a). 

The magnetic current density J* and charge 
density p* can be introduced into Maxwell’s 
equations in the same way as the electric charge 
and current, viz., 


div H=4rp*, (11) 


10H 
curl E+-——= —4nJ*. (11a) 
c at 


These equations are identical with Stratton's 
Egs. I and III, reference 1, p. 464, except for the 
units used: We are using Gaussian, non-rational 
units, and we measure p* in “‘magnetostatic” 
units, J* in ‘‘magnetoelectric” units. The con- 
tinuity equation corresponding to (11), (11a) is 


1 dp* 
div J*+-—=0. (12) 
c at 


We shall assume the time dependence of all 
quantities to be as 
(12a) 


so that (3a) represents an outgoing spherical 
wave. Then (12) becomes 
div J* =ikp* (12b) 
with 
k=w/c. (12c) 
It need hardly be pointed out that J* and * 
have no physical meaning. 

We shall not use magnetic volume currents 
and charges but only magnetic surface currents 
(density K) and surface charges (density 1); 
n corresponds to a discontinuity of H, at the 
surface [cf. (11)] and K to a discontinuity of 
Eten [cf. (11a)]. The quantities K and 7 satisfy 
a continuity equation exactly like (12b), vis., 


div K =ikn. (13) 


The electric and magnetic field can be ex 
pressed in terms of K and 7 most conveniently 
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with the help of a scalar and vector potential, 


vi3., 
E=curl F, (14) 
10F 
H=-—-grad y (14a) 
c ot 


Equations (14) and (14a) automatically satisfy 
the Maxwell equations 


10E 

curl H=- —, (15) 
c Ot 

div E=0. (15a) 


From Eqs. (11), (11a) we find y and F in exact 
analogy to the electric case. If only surface 
charges are present, we obtain 


F()=- f (16) 


vin) = f (16a) 


Inserting these results, and the time dependence 
(12a), into (14a) we find the magnetic and 
electric fields explicitly [cf. Stratton, reference 1, 
p. 466, Eqs. (19) and (23) ] 


= f K(r')Xgrad gdo 
and [Stratton, reference 1, Eqs. (20) and (23) ] 
f grad g)de. (18) 


The gradients are taken with respect to the 
coordinates of the field point r (cf. Section 2 for 
differences between our and Stratton’s notation) ; 
the integral goes over the area of the hole. 

It is seen immediately that (17) satisfied the 
boundary condition (8a), viz., Eran vanishes on 
the screen everywhere outside the hole. In the 
hole, Eran is of course not zero but directly 
related to K [cf. Section 7, Eq. (61) ]. However, 


the problem remains to determine K in such a | 


way as to satisfy the boundary conditions for 
the magnetic field in the hole, viz., (10), (10a), 
and this is, of course, a more difficult task. 
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a. Determination of y from His 


In the solution of our problem, we are greatly 
helped by the fact that the hole is small. Then 
the retardation may be neglected, and H may 
be considered essentially as the magnetostatic 
field corresponding to the charge distribution ». 
The first term in (14a) or (18) is small; indeed, 
from (13) we find that K is of order nka where a 
is the radius of the hole; moreover, grad ¢ is of 
the order ¢g/a so that the first term of (18) is of 
the order of (ka)? times the second. Since the 
hole is assumed to be small compared with the 
wave-length (ka<1), we shall neglect the first 
term of (18) altogether. Then (14a) reduces to 


H= —grad (19) 


and since we also neglect the retardation in ¢, 
(16a) reduces to 


dr’ 
vin = f (20) 


In our approximation we may assume Huan to 
be constant over the hole and equal to 4H, 
(10) ]. Then we have from (19) 


torr. (21) 


The potential problem given by (20), (21) is 
fairly well known from electrostatics ;* we seek a 
two-dimensional charge distribution which gives 
a constant field, Ho, inside the region occupied 
by the charges. It is known that a constant 
inside field is produced by a uniform distribution 
of dipoles in an ellipsoid, the dipoles having (in 
simple cases) the same direction as the field. If 
we now assume the x axis of the ellipsoid 2h to 
be very small, the ellipsoidal charge distribution 
will be equivalent to a surface charge distribu- 
tion. The cross section of the ellipsoid in the YZ 
plane should, of course, be taken equal to the 
hole so that we obtain a rotational ellipsoid with 
semi-axes, a, a, h. The surface density of dipoles 
is proportional to the ordinate of the ellipsoid, 
i.e., proportional to 


p= 


* The following ment wes taken from H. Hertz’s 
solution of the elastic problen of the contact between 
two spheres. See H. Hertz, Crelle’s J. 92 (1881) or A. E. H. 
Love, Theory of Elasticity. 


(21a) 
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The surface charge density is then proportional to 

| Ho-grad tll 22 

o-grad p= (22) 

It would also be possible to determine the 
coefficient of (22) from electrostatics. 

We shall, however, determine this coefficient 

by direct calculation, thus at the same time 


. verifying the solution (22). We put therefore 


(23) 


In the integral (20), we introduce the coordinates 
p=|r—r’| and 8, the angle between the vectors 
r'—r and r (cf. Fig. 2). Then if @ is the angle 
between r and Ho, we have 


Ho-r' =Ho-r+Hop cos (a—£). (23a) 


The integral Eq. (20) becomes, when we insert 
(21) and : 


cf 


cos (a—8)]=}Ho-r. (24) 


The factor p from the element of area cancels the 
factor g=1/|r—r’'| in (20). It is more convenient 


_ to integrate in (24) along a whole chord such as 


RS in Fig. 2, instead of separately over the parts 
RP and PS. Then the integral over 8 goes only 
from 0 to z, and assumes positive (on PS) as 
well as negative (on RP) values. 

Now we find from Fig. 2, 


a@—r?=a (25) 
where s=QS is half the length of the chord and 
x=QT=p+QP=p+rcosB. (25a) 


x goes from —s to +s, and we have 


& 
(26a) 


Then (24) becomes 


Hor cos 8) cos (a—6) 
=3Ho-r. (27) 
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Fic. 2. Illustrating the integration, (24) t 
O is the center of the hole, r=0, P is th r field 06, 
RS any chord through P, T any point r’ on this chord, 


Integration over 8 gives 
cos a)=3Cr*Hy-r (27a) 


so that 
C=1/r*. (27b) 


The assumption (23) therefore actually solves 
the integral equation, and the magnetic charge 
density is 
1 Hey’ 
ovr’. (28) 


From the charge density we can obtain the 
current density K, using (13). We find 


tk 
K =—(a?—r'*) (29) 


Remembering that Ho has vanishing normal com- 
ponent everywhere on the plane x=0, we see 
that K is entirely in the plane of the hole, as it 


should be. 
b. Boundary Condition for £,, 


The correct normal component of Eo, wis, 
E2,=}Eon, must also be obtained by a suitable 
distribution of magnetic currents K, from (17). 
It can easily be seen that the distribution (29) 
does not contribute appreciably to (17). Expres 
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sion (29) is of the order kaHo, grad ¢ is of order 
{/a?, and the integral (17) goes over an area of 
order a?. Therefore the contribution of (29) to 
(17) is at most of order kaH while actually the 
normal component of Ep is of the same order 
as Ho. [In reality, the spatial dependence intro- 
duced by the grad operator in (17) makes the 
contribution of (29) of the order (ka)*Hp only.] 
Therefore we must obtain an additional mag- 
netic current distribution Kg to fit the boundary 
condition for E,. However, we must take care 
that the new current distribution does not give 
rise to an additional magnetic charge density 7 
which might destroy the agreement obtained 
above for the magnetic field. This is most easily 
achieved by letting the magnetic current lines be 
closed circles. Mathematically we must have 


div Kre=0. 


For the actual evaluation it is most convenient 
to use (14), (16). The given electric field E is 
normal to the surface and had the value }£o,, 
which may be considered constant over the hole. 
Therefore (14) gives 


F=}EoXr. (30) 


The integral equation for the components of K 
is the same as for 7 in Section 2, e.g., 


dy'dz’ 
f (31) 


which corresponds to (20), (21). Therefore we 
have, in analogy to (28), 


1 
K,= Eo, 22 


(31a) 


and altogether for the current giving the required 
value of E. 


1 
Kr= r' XEo. (32) 


It is easily seen that this solution satisfies the 
condition (30a). Moreover, since Ey is normal to 
the plane, the vector Kg is entirely in the plane 
of the hole as it should be. The contribution of 
Kz to the first term in (18) can be shown to be 
of order (ka)?Eo, which is negligible. 


c. Final Formula 


The total magnetic current and charge density 
in the hole is now 


1 
K =KytKe r'?)4H 


(33) 
2(a?—r’?)4 


1 


(380) 


5. CALCULATION OF THE DIFFRACTED FIELD 


The diffracted field E, H for x>0 can be 
calculated by inserting (33), (33a) into (17) and 
(18). According to its construction, this field will 
satisfy Maxwell’s equations and all boundary 
conditions. 

We shall carry out the integration at large 
distance from the hole, kr>1. Let « be a unit 
vector in the direction of r, i.e., in the direction 
of propagation of the diffracted wave. Then we 
have 


grad (34) 


In some of the integrals we may replace ¢ by 
its value for r’ =0, viz., 


(35) 


In other cases, it is necessary to expand ¢ in 
powers of r’ and keep the linear term; then we 
obtain 


¢= go(1—tk«-r’). (35a) 


a. Evaluation of the Electric Field 
According to (17) and (34), the electric field is 


E(r) = —ikex Kedy'ds’ (36) 


It is convenient to calculate separately the con- 
tributions of Ky and Kg [cf. (33) ]. We have 


ik a 
f Knedy'ds! =—Hoge [ 
0 


(37) 

2% 

=—ka*Ho¢go. 
3nr 
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The contribution of Kg is slightly more difficult 
to calculate because Kg contains the factor r’ 
and therefore ¢ must be expanded in the form 
(35a). The first term, go, does not give any con- 
tribution so that we obtain 


tk go 


2r? 


where 7’ and @ are polar coordinates in the plane 
of the hole. Remembering that Ep is in the x 
direction (normal to the hole), we find 


f (37D) 


f 


and 
1 
f K ggdy'dz’ K@o. (38) 
us 
Inserting (37) and (38) into (36) we find 
1 
E= (2H o+EoX x). (39) 


We see that E is always perpendicular to x, 
i.e., the waves are transverse. Moreover, for a 
point r on the conducting screen, « lies in the 
plane of the screen and Hy and xX Ep do likewise. 
Therefore E is normal to the screen thus satis- 
fying the boundary conditions. 


b. Magnetic Field 


The first term of (18) follows directly from 
(37), (38) : 


=ik f Kede 
(40) 


= yo(2Ho+Eo Xx). 
3x 


This field is not transversal; only the second 
term (EoXx) has this property. From (39) and 
(40) we see that 

E=-«XH™., (40a) 


The second term of (18) gives 


Hi =— fy grad — f (41) 
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Inserting (33a) we get 


k? 


Now 
f 
and therefore 
H® “sities Ho-«. (42) 
3x 
Remembering the vector identity, 
Ho—Ho:« x=xX(HoXx), (42a) 


we find the entire magnetic field, (40)+(42), 
becomes 


1 
H=H+H® = oor 
X(2HoX«—Eo). (43) 


This expression is a transverse wave as it should 
be. Thus the field H® serves to satisfy the 
transversality condition. In fact, since (42) is 
entirely longitudinal (in direction of x), H® 
serves only to eliminate the longitudinal component 
of H™ and leaves the transverse component un- 
changed. This result will be important in the 


application to the cavity problem (Section 9). 


Comparing (43) and (39) we see that 
H=xXE (44) 
E=-«XH (44a) 


as is required by Maxwell’s equations for plane 
waves. 


c. Total Radiation 
Poynting vector of the diffracted field is 


kta 
Ho — KX«XEo)*. (45) 


Now let 6 denote the angle between the propagz- 
tion vector x and the vector m normal to the 
screen (x direction) and @ the azimuth of x, ie, 
the angle between the plane of « and m and the 
plane of Hy and n. Then the intensity of radiation 
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in the direction x is, per unit solid angle, 


+(sin 0Eo,.—2Hosin a)?]. (46) 


The radiation is thus not symmetrical about any 
axis. The total radiation in all directions is 


f sin oa f dar*| S| 
0 0 


= k*a*(4Ho?+ Eo’). 
7? 


(47) 


The expressions Ey? and H,* on the right-hand 
side denote the time averages. 


6. DISCUSSION OF THE RESULT 
a. Comparison with the Kirchhoff Solution 


‘The result for the diffracted field is entirely 
different from that of the Kirchhoff theory. 
Quite obviously, the polarization of the radiation 
is different because our solution satisfies the 
boundary condition while Kirchhoff’s does not. 
Also the angular distribution of the total radia- 
tion is different (cf. Subsection c, below). How- 
ever, the most striking difference is the absolute 
value of the field. As is easily seen from expres- 
sion (5), the Kirchhoff solution gives magnetic 
and electric fields of the order 


0- 
Our solution gives instead [cf. (39), (43) ] 


which differs from Hx by a factor ka. In addition, 
our solution has a smaller numerical factor. 
Therefore, for small holes, the radiation trans- 
mitted through the hole is very much smaller 
than Kirchhoff’s theory would indicate. The 
fields E and H are reduced roughly in the ratio 
a/\ where is the wave-length of the radiation. 
The radiation intensity is therefore reduced by a 
factor of the order (a/d)?. 


b. Representation by Electric and 
Magnetic Dipoles 


Turning now to our solution itself, we note 
that the field may be considered as owing to an 
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electric and a magnetic dipole. The field of a 
magnetic dipole of moment M at the origin is, 
according to usual electromagnetic formulas: 


H=k'goxX(MXx), (48) 
(48a) 


Similarly, the field of an electric dipole of mo- 
ment P is 


E=kgoxX(P Xx) (49) 
H=kgox XP. (49a) 


Comparison with (39) or (43) shows that the 
hole is equivalent to a magnetic dipole 


2 
0 (50) 
3a 
and an electric dipole 
1 
P=——a'E,. (50a) 
3x 


Each of the dipoles is antiparallel to the re- 
spective field in the hole as might be expected. 
Accordingly, the magnetic dipole is in the plane 
of the hole and the electric one in the direction 
of the normal. These directions are just those 
required to fulfill the boundary conditions for the 
electric field on the screen. The magnetic moment 
(50) can easily be calculated from the expression 
(33a) for the magnetic charge density; we have 


f nr'dr’ 
fr. 
—_ Hat. (51 


It will be remembered that the solution in 
Section 4a was actually obtained by assuming a 
distribution of magnetic dipoles in a flat ellipsoid. 
As is well known from electrostatics, the inner 
field produced by such a distribution is propor- 
tional to the polarization per unit volume (for 
any given shape of the ellipsoid); the total 
moment is therefore proportional to Ho times the 
volume which in turn is proportional to a’. 
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Dimensional considerations, in combination with (m2 and m3 integers), we have Hy,=H),=0 and Th 
the fact that the relation between Mand Hy isa Then the radiated intensity [cf. (46)] is thi 
purely magnetostatic problem, also show that proportional to sin® @ where 0 is the angle between 
the moment M must be proportional to Hoa’. the normal to the screen (direction of Ey) and 

The different factor (1 instead of 2) in the the direction of propagation. This intensity is 
electric moment (50a) is no doubt owing to the maximum on the screen and zero in the normal 
fact that the electric moment is in the direction direction, quite different from the previous case. 
of the small axis of the ellipsoid while the mag- 
netic moment is in the direction of one of the c. Diffraction of a Plane Wave 
We coder of magnitude of M and P isin . it is of some interest to discuss the case of an 
general the same although there are cases for incident plane wave. Let Hi, Es, « be, respec. a 
which one of them vanishes. This is frequently tively, the Sold, the electric field and 
the case for the electric moment P: It is only Ange ‘ brand in the direction of propagation for Rox 

he incident wave. Analogously to (44), (44a 

necessary to assume incident plane waves with ' wr ) : 

we have sin 
the electric field perpendicular to the plane of Al 
incidence. Then the electric field Eo will have no H;=«iXE,, (53) a 

‘component perpendicular to the screen, and only E H ase 
Ho will be different from zero. The radiation (53a) 
intensity is then proportional to sin? x where x is_ and the Poynting vector is 
the angle between the direction of Hp» and the 
direction of propagation of the diffracted waves x. c c 
This radiation is large in the normal direction, seit. salted (54) | and 
and also on the screen in the direction per- (2x) 
. pendicular to Ho, i.e., in the direction of the We shall denote by # the angle of incidence, | 
electric vector of the incident wave. ie., between x; and m (the vector n is in the x 
To make Hp vanish and E)+0, it is necessary direction); then 0<0<x/2. We have to distin- 
to assume standing waves on the “left-hand  guish two cases. 
side”’ of the screen. For example, in a rectangular Inse 
cavity with the sides L,L2L3, the relevant com- a. Electric Field of Incident Wave Perpendicular cros 
ponents of the electric and magnetic fields are to Plane of Incidence P 
given by For this case we have E; perpendicular to the a 
plane of «; and Then has no normal com- an 6 
E,=A, cos ; sin : sin (52) ponent, and the tangential component of H; is . 
1 2 3 
Hi; tan=H; cos 8. (55) 
H,=B, cos sin cos——, (52a) This tangential component is doubled upon re- 
Ly Ly Ls flection so that Hy=2H;; tan; it lies in the plane 
aiid ama ™ of incidence. If, then, 8 is the angle between the and 
ae ees emer Gp sin , (52b) Plane of « and m and the plane of x; and n, we 
Ii Le Ls have for the diffracted radiation intensity per 
unit solid angle 
where A,, B,, and B, are constants and my, m2, m3 Pm 
mee = 5 cos* 6 cos* A) Pie 
m 
and (52c) * =—k‘a®S; cos? #(1—sin? @cos? 8). (56) for 
Ls Or? the 


56) 
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The “diffraction cross section” of the hole for 
this polarization is then 


fisie sin 


= =—k‘a* cos? 3. (57) 
A S; 27x 


B. Electric Field in Plane of Incidence 


If we choose the plane of incidence as X Y 
plane, then ; has the components (cos #, sin #, 0) 
and E; has the components E,(—sin 3, +cos #, 0). 
According to (53), H; has then the components 
E,0,0,1). The components of « are (cos 8, 
sin 6 cos 8, sin @ sin 8). Taking into account the 
reflection and the fact that only the normal 
component of E£; is effective, we get 


—2E,(sin 8, 0, 0), 
Hy=2E,(0, 0, 1), (S7a) 


and [cf. (45) ] 


(2xX Ho— «X« X Eo)? 
=4(xX Ho)? +(«X Eo)? +4x- (Ho X Eo) 
=4E?2[4(1—sin?® 6 sin? 8) +sin? @ sin? 3 
—4sin @sin d cos 8). (57b) 


Inserting into (45) we find for the differential 
cross section 


dA S| 
sin ededB 


16 
=—k‘a*[cos? 0(cos* B 


+i sin? 3)—sin @sindcosB}] (58) 


and for the total cross section 
64 
=——k‘a°(1+} sin? 9). (59) 


The most remarkable feature of the angular 
distribution (58) is the dependence on the azi- 
muth 8: For a given 6, the radiation is smaller 
for 8=0 (i.e., in the direction closest to that of 
the incident wave), than in the “reflected” 


direction 8 =x. We suspected an error in sign to 
be the cause of this result, but we have found 
none. Perhaps the most convincing check on the 
sign is provided by Eqs. (48) to (50a). The result 
mentioned is of course exactly the opposite of 
that expected from any elementary considera- 
tions based on the Huygens principle. For any 
value of 6, the radiation is a minimum in the 
azimuth given by 


cos By =} sin 3/sin 6. (58a) 


In the plane of the screen, the radiation is zero 
for 8 =Bo, which in this case lies between 60° and 
90°. 
The total cross section (59) is greater than 
that for the a-direction of polarization, (57), 
except for normal incidence when the two ex- 
pressions are equal. For unpolarized light, 


64 
A’ sin? #). (59a) 


All cross sections are proportional to \~‘ as in 
Rayleigh’s theory of the scattering by small 
objects. Also the proportionality with a® is the 
same as in the scattering theory. The cross 
section is of the same order as the scattering 
cross section of a dielectric sphere or disk of 
radius a and dielectric constant of the order of 2. 


d. Comparison with Sommerfeld’s Solution 


One characteristic feature of our solution is 
that the current and charge distribution K, » de- 
termine the field only in the half-space on the 
right of the screen. The perturbation of the 
field on the left, H, and E, (cf. (7), (8) ] would be 
obtained by a current and charge distribution 
—K, —»n. This reversal of sign is unusual in 
electrodynamic theory. But a somewhat similar 
procedure had to be used by Sommerfeld in the 
wedge problem. In order to take into account 
that the field is entirely different on the two 
sides of the screen, Sommerfeld put a ‘“‘branch 
plane”’ in the screen and obtained the solution of 
the problem in a Riemann space. This is the 
mathematical expression for the fact that the 
field is not given by the same integral representa- 
tion on the two sides of the screen. 
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7. FIELD IN AND NEAR THE HOLE 
a. Inside the Hole 


From the conditions (10), (10a), (7), (9a), 
(9c), we see immediately that in the hole 


Han tan; (60) 
= Eo, ne (60a) 


Han and E, are thus halfway between their 
values in the “unperturbed” fields on the right 
and on the left side of the screen. On the other 
hand, those field components which would be 
zero if the hole were absent, i.e., the normal com- 
ponent of H and the tangential component of E, 
are determined directly by » and K. Since 4x7 
measures the discontinuity‘ of H,, we have, at a 
point just to the right of the hole, 


(60b) 


2 

In the hole, therefore, the normal component of 
H is of the same order of magnitude as the 
tangential component. This does not in any way 
invalidate our calculations because we have not 
used a perturbation method in which H, was 
considered small compared with Hiss, but we 
have satisfied the boundary conditions exactly. 

In contrast to Hian which was assumed to be 
constant, H, varies rapidly over the surface of 
the hole. Near the boundary r=a, H, becomes 
infinite as 1/(a—r). An integrable infinity of 
this type occurs frequently at edges. It occurs, 
e.g., in Summerfeld’s solution of the diffraction 
by a half-plane, in the neighborhood of the 
diffracting edge. 

The tangential component of the electric field 
is given by [cf. Stratton, p. 467, Eq. (23) ] 


Evan =20nXK, (61) 
where n is in the positive x direction (Stratton’s 
n is in our negative x direction). Inserting (33) 
and remembering that in the hole r-n=0, we 
have 


2 
(61a) 


(a?—r?)3 
‘The discontinuity of H, at x=0 is, of course, only 
mathematical. As explained in Section 6d, the current 


distribution K, 7 determines the field only for x>0. The 
physical field H, is continuous at x=0, according to (9b). 


The second term is in general negligible com. 
pared with the first since H» and Eo, are usually 
of the same order of magnitude. The first term 
is again of the same order as Eoz [cf. (60b)] and 
is singular near the edge of the hole. The tan: 
gential electric field is directed radially outwards 
from the center of the hole. 


Outside the hole, Ets and H, are of course 
zero. The other field components have singy. 
larities near the edge, behaving as 1/(r?—g?)}, 
At larger distance from the hole, His, and E, 
behave as Hoa*/r*, up to distances of the order }. 
For still larger distances, the solutions of Sec. 
tion 5 hold. The field in the neighborhood of the 
hole has not yet been calculated in detail. It js 
clear, however, that His. and E, are very far 
from zero on the screen in contrast to the 
assumption made in Kirchhoff’s theory. 

The tangential component of H is especially 
important because it determines the dissipation 
of energy (joule heat) in the screen if the latter 
has finite conductivity. The energy dissipation 
per unit time is 


f A, (62) 
where v=w/2z is the frequency, 
1 
2x(cv)! (62a) 


is the thickness of the current carrying sheath in 
the conductor, ¢ is the conductivity in electro- 
magnetic units (=10-® times conductivity in 
mho), and the magnetic permeability has been 
assumed to be 1. The integral in (62) goes over 
the whole surface of the conductor. Since Hin 
behaves as 1/(r?—«a*)! near the edge of the hole, 
(62) diverges logarithmically ; it must be cut off 
at a value of r—a of the order of D. The power 
dissipation on the right-hand surface of the 
screen becomes then of the order 


P~H¢vDa? (a/D). (62b) 


An accurate calculation of this quantity would 
be useful for the theory of cavities with small 


holes. 
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8. POSSIBLE EXTENSIONS OF THE THEORY 


In its present form, the theory is only appli- 
cable to holes small compared with the wave- 
length. However, it seems possible to extend it 
to holes comparable in size with the wave-length. 
In this case, the given fields Hy and Ep will 
contain factors of the type e*-’’. Similarly, the 
variation of the factor e#!*-r’! in the Green's 
function g must be considered. From an approxi- 
mate consideration, it seems that the correction 
terms will be of relative order (ka)? rather 
than ka. 

It seems certainly possible to obtain solutions 
in terms of power series in y’ and 2’. The integral 
(20), for instance, can be evaluated if » is any 
odd power of (a?—r’*)!, giving for Ho, a power 
series in r?. From this type of solution, others 
can be obtained by differentiation with respect 
to y and z. We believe that in this way a solution 
for an arbitrary function Ho,(y, z) can be ob- 
tained. However, a more elegant solution may 
perhaps be found in terms of the electric oscilla- 
tions of an ellipsoid. 
~ A second question is the dependence on the 
shape of the hole. From the way in which our 
solution was obtained in Section 4a, it is clear 
that a solution can also be found for an elliptical 
hole. The solution for an ellipse of arbitrary 
eccentricity should give a sufficient idea about 
the dependence of the diffraction on the shape of 
the hole. However, the case of a rectangle may 
also be solvable. 

A more difficult question is the extension to 
screens and holes which are not plane. The case 
of a small hole in a curved surface will probably 
still give a result very similar to ours, as long as 
the surface may be considered plane over a 
region large compared with the hole. This is 
usually the case since the radius of curvature of 
cavities is ordinarily of the order of the wave- 
lengths. If this condition is fulfilled, the boundary 
condition on the surface will be violated only in 
the wave zone where the diffracted field is no 
longer very large. Perhaps a method of suc- 
cessive approximations will be applicable to this 
case. 

However, a different situation arises when the 
hole itself is essentially curved. A simple example 
is a gap in a cylindrical wave guide; in this case, 


the gap itself is the surface of a short cylinder. 
Such a problem will presumably require an 
entirely new solution, but it is likely that at 
least the same principle will work as in our case, 
and that the symmetry of the cylinder will be 
helpful for obtaining the solution. 


9. APPLICATION TO THE THEORY OF CAVITIES 


Condon‘ has given a most convenient theory 
of the excitation of cavities by electric current 
loops. In this theory, the vector potential A is 
developed in terms of normal modes A,, thus: 


A= Lim Pult)An(r). (63) 


By inserting this exprdssion into Maxwell’s equa- 
tions, a differential equation is obtained for p,, in 
terms of the current J exciting the cavity, viz.: 


Wm dpm 


+ ds, (64) 
dt? Qn dt Wm" Pm V m ’ 


where w, is the frequency of the mth normal 
mode, Q» is the dissipation constant, V is the 
volume of the cavity, and the integral is taken 
along the conductor carrying the current J. The 
current is measured in electromagnetic units. 

Since our solution is obtained in -terms of 
magnetic currents and charges, it is convenient 
to expand the field in a slightly different way, 
by using a magnetic vector potential F as defined 
in (14), (14a). The scalar potential y will be 
unimportant because the magnetic field is trans- 
verse. We expand: 


(68) 
and have 
dt 
E= Le QmkmA m(1) (65b) 
with 
Aw=curl Fa/Rm- (65c) 


The F,, form an orthonormal set: 
f F,- Fad V = Voim, (65d) 


where V is the volume of the cavity. Our g,, and 


5 E. V. Condon, J. App. Phys. 12, 129 (1941). 
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Condon’s pm stand in the relation of coordinate 


and momentum. 
The solution (65) satisfies automatically the 
first Maxwell equation (15), (15a) (cf. Section 4). 
The F,,’s are determined in such a way that 


also 


div F,, =0. (66) 
Moreover, the F,, satisfy the wave equation 
with 

Rn =Wm/C. (66b) 


Because of (15), (66), (66a) the second Maxwell 


equation 
10H 


curl E= —-— (66c) 


c Ot 


is also fulfilled by the particular solutions F,,. 
Now assume a distribution of magnetic cur- 
rents of current density J* per unit volume. Then 
the Maxwell equation (66c) is replaced by (11a). 
Inserting Eqs. (14), (14a) we obtain 
1 
curl curl F+— ——= —4nJ*. (66d) 
at? 


Using the expansion (65), the condition (66), and 
the wave equation (66a), and expanding the 
right-hand side of (66d) in terms of eigen- 
functions F,,, we obtain 


4xc? 
fdVI* Fe (67) 
dt? Q,» dt V 


The second term has been added to take into 
account the energy dissipation in the walls of the 
cavity. If we have magnetic surface currents K 
instead of volume currents J*, the right-hand 
side of (67) is replaced by 


dcK- Fa (67a) 


Equations (67), (67a) are obviously the exact 
analog of Condon’s Eq. (64). 

It will be noted that (67a) depends only on K, 
not on ». Indeed, as we have shown in Section 5b, 
the magnetic charge 7 serves only to eliminate the 
longitudinal component of the magnetic field and 
does not affect the transverse component at all. 
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Now 


the normal modes of the cavity all have 


transverse magnetic fields. Therefore to deter. 
mine their excitation we need only the magnetic 
current distribution K, not 7. Moreover, our ex. 
pansion in transverse modes will give us the 
complete field because the longitudinal field pro. 
duced by K is canceled by the effect of 7. 

We may now evaluate (87a) in terms of 77) ,,, 
and Eo,. Inserting (33), we find 


ik 
f Ku- Fado =—Ho- Fa(0) f 


=—ka*H,: F,,(0), (68) 


where F,,(0) is the value of the “eigenfunction” 
F,, at the center of the hole. To evaluate the 
contribution of Kz, Eq. (33), we must expand F, 
in powers of r’ and keep the linear term: 


my my a ° (68a) 


The constant term in (68a) gives no contribution, 
and we obtain 


Jo 0 


OF my OF my 


OF mz OF mz 
ay 


dz’ 


1 
= curl, Fn 


1 
= 2A mz(0). (69) 


In the last transformation, we have used relation 
(65c). Inserting (68), (69) into (67) we obtain 


de 


@m 


— 


Qm dt 


4 
E me(0) —2ikHo-Fn(0))}. (70) 
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We must now determine JJ) and Eo. The 
simplest case, perhaps, is that of a cavity com- 
municating with outside space through a small 
hole. If radiation comes from the outside and if 
there is no radiation in the cavity, then Ho, Eyare 
simply the incident magnetic and electric field at 
the hole. If there is also a field in the cavity, it is 
easy to see that in order to satisfy the boundary 
conditions in the hole, we must choose Hy and Eo 
equal to the difference between the outside and 
the inside field. © 

A case of some interest is that of a cavity with 
a field inside, radiating into space through a 
small hole. Then Ho, Eo must be replaced by the 
inside field at the hole, with negative sign. Some 
caution is necessary since it has been assumed in 
the derivation in Section 4 that the fields depend 
on time as e~*“*. Therefore, if the inside field H is 
to be obtained from (65a), we must put® 


(1/¢)(dqm/dt) = —tkgm. (70a) 


Now let us suppose that the frequency is close to 
the natural frequency w; of a certain mode, 1, of 
the cavity. Then we may put in sufficient 
approximation 


Eo: = —qikiA Izy (71) 
1 
Hy= (71a) 
c dt 
and obtain 
Wm dm 
Q, dt 


4a 
—wigil2wF (0) F,, (0) 
3 V 


WmA 12(0)A mz(0) J, (72) 


where the argument 0 denotes the value of the 
respective function at the center of the hole. 
Equation (72) describes the excitation of all 
normal modes m by one normal mode, 1, through 
the action of the hole. In the stationary state, 
i.e., when 
q1 =A sin wil (73) 


* It can be easily verified that the correct final result is 
obtained by this method even if gm~sin wt or cos wt. It is 
only necessary to go through the derivations from the 
beginning with a different time dependence. 


the excitation of the mth mode will be given by 


4 a @1 
Qm =—-A— sin wit 
V Wm? — w1 


X [2w1Fi(0) Fm(0) (73a) 


Thus quite similar to the perturbation theory of 
wave mechanics, the other modes m will be 
‘mixed in’ with the main mode on account of the 
“perturbation” by the small hole. However, it is 


in general not necessary to use the coefficients gm 


because the field distribution near the hole is 
given more directly by the integrals (17), (18) 
and is only very little modified by the presence of 
the rest of the cavity’ (for given Ho, Eo). 

The most interesting result is obtained by 
applying Eq. (72) to the main mode m=1 which 


gives 


wi dq 


where the emission coefficient y is 
2a’ 2 
= —A;,(0)]. (74a) 


The right-hand side of (74) will cause a shift of 
frequency from w; to 


wi’ =w:(1—7). (75) 


The amount and sign of this shift depends on the 
value of F, and A, at the hole. If A;, is small, the 
frequency will be decreased ; if Fi<A 12, there will 
be an increase of frequency. 

Both alternatives may occur in practice: E.g., 
in a rectangular cavity with electric field in the 
x direction, the field will be given by (52) with 
m,=0. Then the electric field will be zero on the 
XY and the XZ plane; a hole in either of these 
planes will therefore lower the natural frequency. 
On the YZ plane, at the points determined by 


.(52c), Fi:=0 and A,, is a maximum; therefore a 


hole at these positions will increase the frequency. 
If m2=ms3=1, this is the case in the middle of the 
YZ plane (y=43LZ2, z=4L;). For the same mode, 
a hole at y=}L2, s=}L; will leave the frequency 
unchanged. 

At least the first of these results can be 
understood qualitatively from the pattern of 


eter- 
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electric field lines. If the hole is in a plane parallel 
to the undisturbed electric field (XY or XZ 
plane), then in the hole the boundary condition 
E,=0 is relaxed. This is equivalent to a slight ex- 
tension of the cavity, and therefore the charac- 
teristic frequency will be lowered. The electric 
field near the hole will behave about as indicated 
‘in Fig. 3 (exaggerated). 

The case of a hole in the YZ plane is not quite 
so clear. A possible interpretation is that the hole 
reduces the value of E, to one-half and that 
therefore the cavity behaves as if it were smaller, 
giving a higher frequency. However, a closer in- 
vestigation of the fields near the hole seems 
necessary to understand the effects properly.. 

It will be noted that y is purely real. This 
means that only the frequency is changed but not 
the energy dissipation. This is in agreement with 


the result in Section 5c that the intensity of the. 


emitted radiation is of order a® rather than a*. 
The radiation would therefore contribute a dissi- 
pative term of order 7? only. However, it must be 
remembered that our theory does not properly 
take into account the finite conductivity of the 
cavity walls. The dissipation in the walls will be 
increased by the presence of the hole (cf. end of 


Section 7), the increase being of order a?~~7! 


[Eq. (62b)]. 
10. COUPLED CAVITIES 


Let us assume two cavities with one common 
wall with a hole in it (Fig. 1). The fields in the 


| 


Fic. 3. Electric field lines in a cavity with a small hole 
in one of the sides parallel to the field. The density of 
oportional to the field strength, the 


is exaggerated. 


the lines is not 
bulge near the h 


two cavities will be denoted by subscripts a and 
B, respectively. E.g., gam will be the amplitude of 
the mth normal mode of cavity a. Cavity a shal] 
be on the “left” of the boundary x=0, while g 
extends to positive x. We shall assume that mode 
1 of a, and mode 1 of 8 have nearly equal fre. 
quencies, and that the cavities are excited with a 
frequency 


Considering the excitation of cavity 8, we have to 
insert for Hy and Eo, the difference between the 
fields in cavities a and 8. We need only take into 


_account the first mode of each cavity so that, 


similar to (71): 
Eoz=kaigaiA aiz — giz, (77) 
(77a) 


We may drop the subscript 1 and refer in the 
following only to the first mode. Moreover, be- 
cause of (76), we may put kai=kg: in Eo: and Hy. 
Then we obtain the differential equations 


a Ja = —(CaaJa— € 
dt? Q, dt Va 
(78) 
dt Qs dt apd 
with 
Can = = 2F — Ap. 
(78a) 


Cap = Cpa= 2F Fg—AazAsz. 


The values of F,, etc., have to be taken in the 
hole. In all small terms in (78), wa and wg have 
been replaced by wo. Equation (78) has exactly 
the same form as for two conventional coupled 
circuits and can be treated in the same way. 

Let us consider the case of two nearly equal 
cavities, and let us choose the eigenfunctions so 
that the normalized magnetic fields F, and Fs are 
identical at the common boundary. Then, ac- 
cording to (65c) we have also Aagz=Ag-z and 
therefore 

Caa = Cpp = (79) 

We put 


4a’ 2 | 
Y 3 pF (0) A (0)) ( ) 


(78) 
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then y is a dimensionless quantity which meas- 
ures the coupling between the two cavities. Then 
(78) becomes 


ada dqa 


wo dgg 
wo 2 


(80) 


In the last term on the left hand we have kept we 
and wg to allow for a possible lack of tuning. In 
the other (correction) terms we have put 
@a=ws=wo. Neglecting the damping, we find the 
frequencies of the coupled cavities are given by 


w= 4 (wa? + wp") —wory 
(wa? — wg”)? (81) 


If the two cavities are considerably out of tune, 
i.e., if 

| Wa — we | Pwy, (81a) 
they will oscillate almost independently, i.e., 
there will be one proper frequency close to we at 
which mostly the first cavity is excited while the 
second has very little excitation, another proper 
frequency close to ws for which the reverse is the 
case (cf. Section 11). 

The cavities are well tuned if 


| wa — we | (81b) 


The tuning requirement is thus determined by 
the size of the hole, cf. (79). If w2=wg=wo, (81) 
gives 

w? =wo?(1—y+y). (82) 
Thus we get one mode whose frequency is + the 
same as for the uncoupled cavities, 


@1=Wo, (82a) 
and another mode of frequency 
w2=wo(1—+) (82b) 


when we neglect terms of order 7’. 
In the first mode we have [cf. (80) ] 


Ja = Qp- (83) 


This means that the tangential magnetic field at 
the common boundary is continuous, H.=Hs,, 
and the same holds for the normal component of 
the electric'field. The hole, then, does not perturb 
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a. 


Fic. 4. Field distribution in the first mode w; of two 
coupled tuned cavities. (a) Cross section of cavities. 


(b) Plectric field E,, as a function of x (same curve for 
any value of y). (c) Magnetic field H,. 


the field at all and the two cavities oscillate as if 
there were no coupling. This explains immedi- 
ately why the frequency is not disturbed. Never- 
theless, the coupling is important because it 
determines the phase relations between the two 
cavities. For the simple case of two rectangular 
cavities oscillating in their lowest mode we have 
plotted the electric and magnetic field in Fig. 4. 
The electric field is assumed to be parallel to the 
common boundary; it is zero inside the hole just 
as on the conducting wall. 
The second mode has 


qa= — (83a) 


Then, the magnetic field changes its sign when 
the boundary between the cavities is crossed (cf. 
Fig. 5). In the hole, the tangential component of 
the magnetic field is zero. Around the hole, there 
is a large perturbation of the field. The tangential 
electric field has a maximum in the hole. In Fig. 5 
we have plotted the electric and magnetic field 
for the same cavity as Fig. 4, but this time for the 
second mode, w2=wo(1—-y). Since in our case the 
unperturbed electric field has no component 
normal to the common boundary, y is positive, 
and the frequency is reduced by the coupling. 
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zero on the common boundary (x=0) is H,. If we 
eae eee y a. do not restrict our attention to the lowest mode, 
-_ we have for the normalized field : 


32 


b. 


P where m, and ms are positive integers not zero 
Cc. and are the lengths of the three edges of 
y single cavity. Then we have for the coupling 

coefficient [cf. (79) ] 


8 a’ 


Ey ot x- 4 sin? —. (84a) 


(m;/L1)?+(ms/Ls)? Ls 
= The last factor possesses nodes on the common 
LZ \ boundary, x=0, for any value of ms; greater than 


oj Ff H s at Y, one. In general, large coupling is desired ; then the 


hole should be put at or near the maximum of , 
If definite phase relations between the oscilla. 


Fic. 5. Field distribution in the second mode w: of the tions of the two cavities are required, the correct 


cavities of Fig. 4. (a) Cross section [identical with 4(a)]. . 
(b) Electric field Ey along a line not going through the frequency must be selected. It seems possible to 
) Magnetic field along the same line. predict, in many cases at least, whether the 


hole (y= (c 
(d) Electric field E, along a line y= oing through the . ‘i * 
hole. (e) Magnetic field the desired phase relation corresponds to the higher 
or the lower frequency. In any case, two coupled 


This corresponds to the fact that now, so to a 
speak, the two cavities oscillate as a whole; the : 
electric field E, no longer goes to zero in the hole 
but has a secondary maximum there, thus be- 
ginning to develop towards the field correspond- 
ing to the lowest mode in the larger cavity a+8 
with the central boundary absent (Fig. 6). This 

mode is entirely different from that of the sepa- b. 
rate cavities (Fig. 4) ; therefore in mode 2 (which 

is the lowest harmonic for the larger cavity of , a. 


Fig. 6), a large perturbation field exists around 


the hole which varies rapidly in space. The Cc. 
energy dissipation is likely to be considerably rail 
greater than for mode 1 (cf. Section 7, end). 

The coupling coefficient y depends sensitively cil 


on the position of the hole, and may be zero for 


certain positions. For example, in our case of two Fic. 6. (a) The cavity of Fi m4 and 4 without the = 
iti i ; ; in the middle. (b) Electric field distribution in 
rectangular cavities with an electric field in the of tke 6. field in the 


y direction, the only field component which is not — same mode. 
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cavities which are exactly tuned will give two 
modes with a relative frequency difference of the 
order y~a*/V, the phase relations between the 
fields in the two cavities being opposite for the 
two modes. 


11. STEP-UP COUPLING 


It is often desirable to couple two circuits in 
such a way that the field in 8 is stronger than in a. 
This can be achieved for two exactly tuned 
cavities if they are coupled asymmetrically, viz., 
by a hole which is in a position corresponding to a 
high value of the normalized field F, and to a low 
one of Fs (Fig. 7). Then the excitation may be 
expected to adjust itself so that the magnetic 
field is continuous in the hole, viz.: 


Ha=QaFa (hole)=Hs=qeFe (hole), (85) 


where F., Fs are the normalizéd fields. Then, 
since Fs is made small compared with F,, we get 
qa as desired. 

We shall show in this section that (85) is 
- actually fulfilled for one of the two modes, viz., 
that having the same frequency as the uncoupled 
cavities. There is, however, another mode in 
which the field in 8 is smaller than in a. Further- 
"more, we are going to show that two slightly mis- 
tuned cavities coupled symmetrically will give 
the same results as two exactly tuned cavities 
coupled asymmetrically.’ 

The problem of cavities coupled asymmetrically 
in the way indicated in Fig. 7 can be easily treated 
by Eqs. (78) to (78c). It is convenient to define 
the ‘‘displaced frequencies.” 


Wa’? = — (4/3) wo*Caa(a*/ V), 


(86) 
wp’? = wa — (4/3) wocaa(a*/V), 
and the coupling coefficient 
¥ = (4/3) wo*Cas(a*/ V), (86a) 


the volumes of the two cavities having been 
assumed equal. Then, if w denotes the actual fre- 
quency, (78) becomes 


(wa’? — w*) ga = 0, 
— w*)gs=0, 


(87) 
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Fic. 7. Asymmetrical coupling of the two tuned cavities. 


which has the two solutions 
2= 

(87a) 
The ratio of the excitation amplitudes ga, is: 


qa 


(88) 


the upper sign corresponding to the higher fre- 
quency w:, the lower sign to we. d is given by 
/2 w — Wa! 
B 


2ywo? 


(88a) 


Yo 


We shall now make the simplifying assump- 
tions [cf. (78a), (78b) ] that A, is zero in the hole, 
and that F, and Fs have the same direction. 
These assumptions are fulfilled by the cavities 
treated in Section 10 where A = A, and Fis in the 
z direction everywhere on the common boundary 
of the cavities. Then we may write [cf. (78a), 
(78b) } 


Cou ™ Catt, (89) 
Cap =Cap/, (89a) 
with 
u=F,/ Fs. (89b) 
Then (86) reduced to 
Wa’? = Wa" — Yworu, 
(90) 


we’? = wg” — 
and (88a) becomes 
A= (u?—1)/2p. 
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If we now assume that the cavities were originally 
tuned, i.e., wz=ws=wo, we obtain, by inserting 
(90), (90a) into (87a) and (88) 


wi? =wo’, 
ws? 1 —y(u?+1)/u], 
(92) 


(46/qa)2= —1/p. (92a) 


We find thus one mode w; of frequency equal to 
the unperturbed frequencies, and another mode 
we of lower frequency. In the first mode the ratio 
of the amplitudes of excitation gs/qa is such that 
the magnetic fields in the hole are matched, in 
agreement with Eq. (85). This mode, then, gives 
the desired stepping up of the amplitude from a 
to 8. However, in the second mode, the ratio of 
the g’s is reversed (in addition, the phase relation 
is opposite). Caution is therefore necessary to 
excite the coupled cavities with the correct fre- 


(91) 
(91a) 
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quency in order to obtain the desired results, 
(For frequencies other than w: and w2, beats wil] 
occur.) 
Moreover, we see that it is not at all n 
to couple the cavities unsymmetrically in order 
to step up the excitation from a to 8. It is only 
important (cf. (88) ] that A be different from zero, 
ie. [ef. (88a) ], that wa’ and we’ be different. This 
can be achieved either by making wa=wg but 
Caa*Cgg; then we have two tuned but unsym. 
metrically coupled cavities as treated in (89) to 
(92a). But it can equally well be achieved by two 
slightly mistuned and symmetrically (or also 
unsymmetrically) coupled cavities; then waxu 
and Caa=Cgg so that 
@g—-Wa 
(93) 


Yo 


In either case, the more strongly excited cavity ig 
the one whose w’ is closer to the actual frequency, 
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Using Bhabha’s method, the intensity of slow mesons is determined. The intensity of mesons 
of energy greater than 108 ev and less than 2.7 X 10* ev is 6 percent of the total meson intensity. 


Near the geomagnetic equator, the total intensity increases by 54.4 percent, the meson in- 
tensity by 21.7 percent, and the electron intensity by 162.1 percent when the altitude increases 
from 3100 ft. to 7200 ft. The electron intensity, which is about 25 percent of the meson intensity 


at 3100 ft., increases to 50 percent of the meson intensity at 7200 ft. 


HE behavior of cosmic-ray electrons is ex- 

plained satisfactorily by the cascade theory. 
On the basis of this theory, very few, if any, of the 
primary electrons incident on the top of the 
atmosphere can reach sea level. Experimental 
results, however, indicate clearly that a signifi- 
cant part of cosmic radiation at sea level, and 
even at great depths, consists of electrons. It is 
believed that most of these electrons in the lower 
atmosphere arise out of meson decay.’ The test 
of the correctness of any theory of meson decay 
ultimately depends on an accurate determination 
of the proportions of electrons and mesons in 
cosmic radiation. 

One method of determining this ratio of elec- 
tron to meson component is by a statistical ex- 
amination of the tracks obtained in a Wilson 
chamber. For high accuracy, a large number of 
photographs are required. It is not possible to 
distinguish mesons of all energies from electrons 
by their tracks in a Wilson chamber.** Stray 
photographs of low energy mesons may be ob- 
tained but this fact can only be utilized to indi- 
cate the presence of very low energy mesons in 
cosmic radiation. 

The second method is by the use of Geiger 
counters. With a counter telescope, the total 
electron and meson intensity can be determined. 
(The intensity of protons, etc., if present can be 
considered to be negligible and, for purposes of 
this paper, will be included in the meson com- 
ponent.) By introducing an absorber like lead 
between the counters of a counter telescope, 


, asa and Heisenberg, Ergeb. d. exact. Naturwiss 17, 
*Blackett, Proc. Roy. Soc. A165, 11 (1938). 
*S. H. Anderson oe C. D. Nedermeyer, Phys. Rev. 51, 


884 (1937); 54, 88 (1938). 
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electrons can be absorbed and the meson in- 
tensity determined. The difference between the 
total and the meson intensity gives the electron 
intensity. This method has its obvious handicaps. 
If small thicknesses of lead are used to absorb 
electrons, some of the high energy electrons are 
not cut out and are included in the meson com- 
ponent. If large thicknesses of lead are used for 
the purpose, low energy mesons are absorbed and 
they will be included in the electron component. 
Mesons of energy, 3X10* ev can, according to 
Rossi and ‘Greisen,‘ penetrate 21.76 cm of lead. 
Lead absorbers of thickness up to this value have 
been used in estimating the meson intensity. 
Therefore, it is usual to call all mesons of energy 
less than 3X 10* ev as “‘slow’’ mesons. Mesons of 
higher, energy are called ‘‘fast’’ mesons. 

If the absorption curve of slow mesons is 
determined, it is possible to obtain the ratio of 
electron to meson intensity from the results of the 
usual experiments. Besides, the intensity of slow 
mesons of different energies can be determined. 
In recent years, attempts have been made to 
determine the intensity of slow mesons. Auger® 
and Greisen*® make use of the differential absorp- 
tion of electrons and mesons in the same equiva- 
lent thickness of different materials like lead, 
carbon, aluminum, etc., to arrive at an estimate 
of slow meson intensity. Auger’s results are sub- 
ject to considerable error. Greisen’s analysis is a 
very elegant one, and his results are in good 
agreement with the results reported in this paper. 
But such a method is an indirect one, and an 


' bi. 7 Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
1 

*’P. Auger, Phys. Rev. 61, 684 (1942). 

* K. Greisen, Phys. Rev. 63, 323 (1943). 
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ticles. This fact has also been experimentally 
verified by one of us."° 

Any attempt to estimate the intensity of 
‘‘slow”’ mesons must be free from the objections 
indicated. In addition, the results must be free 
from the effects of side showers. In order to 
minimize the possibility of chance coincidences, 
it is desirable to avoid a double coincidence 
system. Impelled mainly by these considerations, 
an’ attempt has been made here to give a report 
of a new approach to the problem made possible 
by the method developed by Bhabha.°® 


EXPERIMENTAL ARRANGEMENTS 


A counter telescope designed on the basis of 
Bhabha’s method is shown in Fig. 1. Counters J, 
2, and 3 were connected in coincidence and form 
the vertical telescope. Counters 4 and 5 were con- 
nected in parallel and in anticoincidence with 
counters J, 2, and 3. The pulse was taken from 
the copper cylinders of counters 4 and 5 for anti- 
coincidence and these counters were shielded. 
Each of the counters used was 6 in. long and 1} in. 
in diameter. This type of counter telescope gives 
results free from the effect of side showers.'° 

Positions I, II, and III were for placing lead. 
High energy electrons produce showers in 1.25 cm 
of lead placed in position II] and are cut out by 
the shower particles tripping the counters 4 or 5 
or both. Lead in position II loses a good deal of 
its significance when lead is placed in position I. 
in. : Even under such conditions, lead was placed in 

position II, for purposes of maintaining identical 

geometrical conditions for all measurements. 
4. Lead placed in position was for absorbing low 
energy electrons. Lead completely covering the 
solid angle was placed in position | for absorption 
measurements. As penetrating non-ionizing cos- 
mic rays form a negligible fraction of the total 


estimate of the slow meson intensity by a more 
direct method is, therefore, desirable. Schein and 
others’? and Sarabhai® cut out electrons while 
measuring the meson intensity by not recording Tanuz I. Measurements of intensity. 
particles which produce showers in small thick- 
nesses of lead. It has been pointed out by Bhabha® ) Time 

a Lea otal in Counts/ 
that a good fraction of low energy electrons pro- Counts hours hour 
duce showers of very few particles, often only of 118412 
one and that this fraction cannot be cut out by 
just not registering the shower-producing par- 


7 Schein, Jesse, and Wollan, Phys. Rev. 58, 1027 (1940). —————— 
8 V. Sarabhai, Proc. Ind. Acad. Sci. 19, 37 (1944). 10Chandrashekhar Aiya, Proc. Ind. Acad. Sci. (in 


®H. J. Bhabha, Proc. Ind. Acad. Sci. 19, 23 (1944). print). 
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SLOW MESONS IN COSMIC RADIATION 


intensity" and as scattering can be considered 
negligible for purposes of the accuracy aimed at, 
this choice of position I for lead in absorption 
measurements is not objectionable. 

The measurements taken are given in Table I. 
From observations 3, 4, 5, and 6, a graph is 
drawn and is shown in Fig. 2. This gives the 
absorption curve of slow mesons. It is extrapo- 
lated to zero thickness by dashes. 


DISCUSSION OF THE RESULTS 


No. 1 represents the total intensity. If we as- 
sume that No. 3 gives the meson intensity, 
the electron intensity is 23.5+1.9 counts/hour. 
Therefore, at Bangalore (geomagnetic ‘latitude 
=3.3°N and altitude=3100 ft. above mean sea 
level) cosmic rays consist of 80 percent mesons 
and 20 percent electrons. It will be clear from the 
discussion which follows that No. 3 gives the 
meson intensity to an accuracy of 2 percent and, 
therefore, the ratio of electron to meson intensity 
stated above is accurate to within 2 percent. The 
difference between No. 2 and No. 3 represents the 
fraction of electrons cut out by splitting the lead 
as suggested by Bhabha. It is 10.6+6.8 percent 
of the total electron intensity. 

According to Bhabha, the fraction of electrons 
which penetrates 5.25 cm of lead and which does 
not produce adequate showers to trip counters 4 
or 5 or both is at the most 2 percent of the total 
electron intensity. Therefore 0.47 count in 
94.8+0.7 counts/hour represents the electron in- 
tensity registered as mesons in observation No. 3. 
This is within the limits of accuracy of the 
measurement. Consequently it is correct to say 
that No. 3 gives the intensity of mesons which 
penetrate 5.25 cm of lead. Observations 4, 5, and 
6 give the meson intensities to a still better degree 
of accuracy. Therefore, the graph drawn with 
these observations represents the absorption 
curve of slow mesons. By making use of this 
graph and the curves of Rossi and Greisen,* we 
can deduce the intensity of mesons within differ- 
ent energy limits. Thus, it is found that slow 
mesons of energy greater than 10* ev and less 
than 2.7 X10* ev is about 6 percent of the total 
meson intensity. 


ase” and Rochester, Proc. Roy. Soc. A181, 399 
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Fic. 2. Absorption of slow mesons. 


One of us has reported something abnormal in 
the behavior of mesons of energy in the region of 
3X10* ev, and the subject is yet under investiga- 
tion. It is, therefore, not proposed to go into the 
question of mesons of this energy here. 

The curve given by Euler and Heisenberg’ for 
the energy spectrum of mesons shows that 
mesons of energy less than 10* ev constitute 
about 2 percent of the total meson intensity. The 
smallness of this fraction is due to the collision 
loss of mesons which increases rapidly below this 
energy. If our graph were extrapolated to zero 


' thickness, which is clearly not justified for the 


reason just mentioned, we would find that about 
2.7 percent of the total meson intensity consists 
of mesons of energy less than 10* ev. This is 
natural. The meson absorption curve for thick- 
nesses of lead less than 5.25 cm must fall off 
towards the X axis, i.e., must be below the dotted 
line in Fig. 2. Although it is not conclusive, it 
may be pointed out that the point for 5.25 cm of 
lead having a higher statistical accuracy than 
others is below the line and can be taken to indi- 
cate that the curve does fall off as is to be ex- 
pected. The percentage of mesons absorbed in the 
first 5.25 cm of lead, i.e., of energy less than 
0.93 X 108 ev is very small (less than 2 percent of 
the total intensity), and it is difficult to estimate 
this small quantity by obtaining readings of 
meson intensity only within the first 5.25 cm of 
lead as any method used does not entirely cut out 
electrons and the fraction of electrons not cut out 
increases with decreasing thickness of lead. 

It would, however, be possible to obtain an 
upper limit to the meson intensity at 1.25 cm of 
lead by omitting the 4.0 cm of lead absorber. 


#2 Chandrashekhar Aiya, Nature 153, 375 (1944). 
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TABLE II. Comparison of results at 3100 ft. and 7200 ft. 


Intensity at Percent increase 


Intensity at 7200 ft.in corresponding 


3100 ft. in counts/hr. to a rise in 

counts/hr. (Ootaca- altitude of 
Component (Bangalore) mund) 4100 ft. 
Total 118.3 182.7 54.4 
Mesons 94.8 121.1 21.7 
Electrons 23.5 61.6 162.1 


Serabhai* has done this with 2 cm of lead. Since 


the intensity of the soft component is 25 percent 
of the meson intensity at Bangalore and we can 
gauge from the figures given by Bhabha® that 
very roughly about 10 percent of the soft com- 
ponent passes through this thickness without 
producing showers, it follows that the soft 
component would contribute about 23 percent to 
the meson reading so obtained. But it is quite 
clear that such a measurement would have no 
meaning if taken at an altitude of 7200 ft. where 
the soft-component intensity is 50 percent of that 
of the hard component and, at still greater alti- 
tudes the error would be much greater. 

The graph in Fig. 2 holds for the solid angle of 
the counter telescope used by us and indicated in 
Fig. 1, which is drawn to scale. For different solid 
angles, a part of the mesons at least will have 
traveled through different distances in the atmos- 


phere and suffered different degrees of losses of 


’ energy. It is, therefore, not quite impossible that 


the absorption curve may be slightly different for 
different solid angles, but this is a matter for 
investigation. 

The counter telescope used for experiments at 
Ootacamund (geomagnetic latitude=1.7°N and 
altitude = 7200 ft. above mean sea level) by one 
of us (S.V.C.) was quite similar to the one used 
here. The results at the two places can be utilized 
to indicate the variation of the different com- 
ponents with altitude near the geomagnetic 
equator. This is done in Table II. 

The electron intensity which is about 25 
percent. of the meson intensity at 3100 ft. in- 
creases to 50 percent of the meson intensity at 
7200 ft. 
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The integration of the Lorentz force equations to give electron or ion paths has been reduced 
to simple quadratures for systems in which the electric field is zero and the magnetic field is a 
function of one Cartesian or cylindrical coordinate. For several interesting types of magnetic 
field variation the quadrature can be carried through analytically; and even for complicated 
magnetic fields, or such as are known only empirically, the numerical integration can be 
effected without difficulty. From general considerations of the functions involved, it is possible 
to determine the extension and periodicity of the orbits for any set of initial conditions. The 
representation used is also convenient for obtaining information regarding the dispersion and 
focusing characteristics of the trajectories, some of which have unique properties of promise 
for use in specific instrument design such as mass spectrometers, beta-ray spectrographs, etc. 
Schematic designs of such instruments are proposed, and a discussion ig given of their ad- 


vantages and special properties. 


INTRODUCTION 
HE focusing and dispersing properties of a 
uniform magnetic field have been applied 
with great success in recent years to mass spec- 
trographs, mass spectrometers, and §-ray spec- 
trometers, as well as to more special uses. 

Many instruments have been built and de- 
scribed which employ either 180° focusing or the 
focusing properties of sector-shaped magnetic 
fields. All are generically of the same class, of 
which the 360° analyzer is the only one having 
the property of perfect focusing. However, it is 
not practical to attempt to utilize the latter 
because of the physical impossibility of locating 
both a source of ions and a collector at exactly 
the same place. 

In contrast to the wide application of the 
focusing and dispersing properties of a uniform 
magnetic field, little use has been made of the 
focusing properties of non-uniform fields. The 
problem of calculating paths, dispersion, and 
focusing characteristics for a general type of non- 
uniform field is one of great mathematical diffi- 
culty. There are, however, some cases of con- 
siderable physical interest in which the magnetic 
field varies in a fairly simple manner, and for 
which it is possible to calculate exactly the 
paths, dispersion, and focusing properties. 

; 

1 While there is a vast literature on the general optics 
of charged particles moving in magnetic fields, applications 


of the type considered here do not appear to have been 
given specific treatment. 


In all cases where the pole faces of a magnet 
or an electromagnet are symmetrically located 
with respect to a median plane, the magnetic 
field at any point on the plane will be per- 
pendicular to it. For an ion or electron moving 
in this plane the problem of describing its path 
reduces to a two-dimensional one. 

For those cases where the magnetic field in 
the median plane can be expressed as a function 
only of.one Cartesian coordinate, or as a function 
only of the radial distance from the axis of sym- 
metry, for cases of circular symmetry, the equa- 
tions describing the ion or electron paths admit 
of a simple treatment. Even if the dependence of 
the field upon such coordinates is known only 
experimentally, the complete solution may be de- 
termined by simple numerical integration to an 
accuracy equal to that with which the magnetic 
field is known. The general theory will be out- 
lined and discussed below and certain cases 
analyzed in detail. 


GENERAL THEORY 


The Lorentz force equation for a charged 
particle in the presence of a magnetic field H 
and the absence of an electric field, is: 


(1) 


where e=charge of particle in e.s.u.; c= velocity 
of light ; and v= velocity of particle in cm/sec. 
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In addition to Eq. (1), we have the equation: 


(4)?+ (y)?+ (2)?=v?=const., 


since the magnetic field can do no net work on 
the charge. For the case of both Cartesian and 
cylindrical coordinates, we may, without loss of 
generality, consider H to be parallel to the z 
direction, so we then have: 


(z)?+(y)? =v", 
(7)?+ (76)? =0°. 


Here v refers to the resultant velocity in a plane 
perpendicular to the z direction. 
Written out in terms of Cartesian components, 


Eq. (1) becomes: 


(2) 


(3) 


where #7 is the magnitude of the magnetic 
field, or: 


#—r(6)?=— 


(4) 


1d 
— —(r26) = 
dt mc 


r 


in cylindrical coordinates. Equations (3) and (4) 
cannot in general be solved explicitly if H is a 
function of both coordinates, as the equations 
will not be linear. 

If, however, H is a function of only one 
coordinate, the equations are solvable in a 
simple manner, as follows: 

Referring to Eq. (3), we may assume, without 
loss of generality, for this case that H=H(x), 
and define a function f=f(x) as: 


=— | H(x)dx. 5 
f — (5) 


Integrating the second of Eqs. (3), we get: 


é 
| Hdx=of. 6 
y —f uf (6) 
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However, from Eq. (2), we have: 
&=+(v?—y*)'= (7) 
Since dy/dx=y/z, we may combine Egs. (6) 


and (7) to obtain: 


ae 
or: 
fax 


We see thus that when H is a function of only 
one Cartesian coordinate, the equations of the 
trajectory are reduced to the simple differential 
equation (8), which permits immediate integra- 
tion. When the magnetic field is a function only 
of the cylindrical coordinate r, the treatment is 
similar. For this case we define a function g(r) by: 


e 
g(r) =——- fear (10) 
vmc r 
Integrating the second of Eqs. (4), we find 
e 
r6=—— | rHdr=vg. (11) 
rmc 


Using Eqs. (2) and (11), we have: 
+0(1—g’)!, (12) 


and since r6/7=r(d0/dr), we have from Eqs. (11) 
and (12): 


dé 1 g 
—=+- (13) 
dr r (1—g?)! 
or: 
gdr 
6= 14 
+f r(1—g?)! 


We see, therefore, that the treatments for the 
Cartesian case and for the cylindrical case are 
similar, and involve essentially nothing more 
than two quadratures, either or both of which 
may be done numerically if not analytically. 

Equations (5) and (10), which define f and g, 
will each contain a constant of integration. 
These constants may be used to specify the slope 
of the trajectory for any desired value of x orr 
by direct application of Eqs. (8) or (13). The 
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(e) (f) (g) 


Fic. 1. Various types of orbits possible in a linearly varying magnetic field, 
determined by the parameter ¢, defined. by fe. (17). 


other initial condition, namely, the starting 
position in the x, y or 6, r plane is taken care of 
by the constant of integration in Eqs. (9) or (14). 
The fact that the equations defining the 
trajectories involve the radical (1—f?)! for the 
Cartesian case, and the radical (1—g*)? for the 
cylindrical case, has a very simple and significant 
interpretation. Since these radicals must repre- 
sent real magnitudes for any physically possible 
trajectory, we see immediately that the extension 
of the trajectory in the x or r coordinate will be 
given by —1<f<1, or respectively. 
If the equation 


1—f?=0 (15) 


has two real roots, the implication is that orbits 
will be bounded as regards extension in the x 
coordinate. This follows from the fact that dy/dx 
will be infinite for x equal to a real root of Eq. 
(15). If Eq. (15) has only one real root the orbit 
is one that comes in from and returns to plus or 
minus infinite values of x. For the case where 
Eq. (15) has two real roots the orbit will be 
either completely closed, or have an infinite 
extension in the y direction. The latter orbits 
will be of periodic behavior and will, in general, 
consist of an infinite number of loops of cycloidal 
appearance. 

The same general considerations as just dis- 
cussed for the case of Cartesian coordinates also 
hold for the cylindrical system and figures demon- 
strating them will follow in the detailed dis- 
cussion of specific cases. 

In the calculations for the paths in either 
coordinate system the quantity e/vmc enters as 
a fundamental constant of the motion. It is 


convenient to define this quantity as a, i.e.: 
a=e/vmc. 


In the discussion of many cases of physical 
interest it is more expedient to refer to the 
potential of a charged particle as a measure of 
its energy rather than to its velocity. We may, 
therefore, write a as j 


a=- 
mV 


where V is the potential of the particle in volts 
or its energy in electron volts. 


LINEARLY VARYING MAGNETIC FIELD 
The linear magnetic field may be defined by: 
H=bx, (16) 


the origin being placed where the field is zero. 
From Eq. (5) it follows that: ° 


(17) 


abx? 
f= a f = é, 


where @ is the constant integration and a=e/vmc, 
as previously indicated. 
Introducing the change of variables: 


so that g and & are dimensionless, we find Eq. (8) 
takes the form: 


(18) 
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By inspection of the right side of this equation 
it is clear that if it is to represent a real trajectory 
é<1. Moreover, the nature of the real trajectories 
will be determined by the magnitude of 2, as 
follows: 

For 1>é¢>-—1, the trajectories will all cross 
the y axis, with a slope = 2/[1 —2*]. They will be 
symmetrical with respect to y, and the iorf or 
electron will oscillate periodically between fixed 
values of x. If 1>2>0, the trajectory will simu- 
late a sinusoidal oscillation with a monotonic 
increase or decrease in y, as illustrated in Fig. 
1(a). When ¢=0, the trajectories will have an 
inflection when crossing the y axis, as shown in 
Fig. 1(b), so that the change in y is still mono- 
tonic as the path is traversed. For 0>¢> —1, the 
variation in y is no longer monotonic, and the 
trajectory tends to form loops as it oscillates 
between the fixed limits in x [cf. Fig. 1(c)]. At 
¢=—0.46, the loops of consecutive periods on 
either side of the x axis touch, and become inter- 
laced [ef. Fig. 1(d)] as ¢ is still further decreased. 
Complete overlapping of the loops, so as to form 
a single figure-eight trajectory, develops for 
@=—0.6522. For the range —0.6522>¢>—1, 
they are again similar to a network of interlaced 
figure of eight oscillations, which spread out and 
ultimately degenerate into the single pair of 
split loops when ¢= —1, as shown in Fig. 1(e). 

In this whole range of @ the maximum exten- 
sion of the oscillations, which may be obtained 
from the roots of Eq. (15), continually increases 
as @ decreases, according to the equation: 


Emax = (1—2)}. (20) 

- Moreover, from the roots of Eq. (15) the width 
of the oscillations is given by: 

Emax — Emin = (1—2)*— (1+2)}. (21) 

The detailed calculation of the trajectories 

may be readily made by introducing the substi- 

tution: £=(1—2)!cos® and reducing the in- 


tegral of Eq. (19) to standard elliptic integral 
form. It is thus found that: 


1 
—K(k)—2E(#, k)+F(#,k)], (22) 


where E(®, k), F(#, k) are the incomplete elliptic 
integrals corresponding to E, K. It is to be noted 
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that Eq. (21) applies directly only to the quarter. 
cycle of a single complete oscillation passing 
through the origin. Its continuation to both 
preceding and following elements of the trajec. 
tory is easily made by simple shifting of the 
initial point and appropriate changes in sign, 
as dictated by Eq. (19). 
The wave-length, along y, of the oscillations jg 
given by: 
d= 2v2 | 2E(k)—K() |, (23) 


where K and E are the complete elliptic integrals 
of the first and second kind, the modulus 
k=[((1—2)/2]!, and the absolute value is used 
to take care of the change in sign of the enclosed 
term at €= —0.6522. dX is expressed in the same 
units as J. It will be noted that \=49(Zmax). 

In the limiting case ¢=—1, integration of 
Eq. (19) readily gives: 


v2+ 


—(2—#%), £>0, (24) 


from which were plotted the trajectories of 
Fig. 1(e). 

The trajectories of greatest practical interest 
are those for ¢< —1. The general features of these 
can also be deduced by direct inspection of 
Eq. (19). These no longer cross the @ axis, but 
split up into two identically shaped individual 
trajectories, for each value of ¢, placed sym- 
metrically on the two sides of the g axis. Each 
of these is composed of a series of interlaced or 
separated loops extending from 


|Zmax!=(—1—2)', to |#min| =(1—2)#, (25) 


and hence recede from the g axis as é is decreased. 
The overlapping of the loops increases as é de- 
creases [cf. Figs. 1(f) and 1(g) ]. The wave-length 
is given by: 


(1-—)x-z] (26). 


where now the modulus k=[2/(1—2) ]. A plot 
of X and (Zmax—Zmin) VS. Emax is given in Fig. 2. 
It will be seen that both \ and the lateral range 
of the oscillations decrease as Zmax increases—i-., 
as the trajectories are confined to regions of 
greater magnetic field. This is, of course, to be 
expected from elementary considerations. 
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For calculating the exact shape of the trajec- 
tories for this range of 2, the integral of Eq. (19) 
may-again be transformed into standard elliptic 
integrals, by setting Z*=cos 2@—¢ with the result 
that: 

g=(1—2)!LE@, k) —(1—4k*) F(®, k)]. (27) 


As remarked above, orbits of the type illus- 
trated*in}Figs. 1(f) and 1(g) are of greatest 
interest. This is because of the practical difficulty 
of obtaining a field variation in which H is an 
odd linear function of x. Also, for any restricted 
region of an actual non-uniform field H may be 


— 


| 2 3 4 5 6 
Xmox. 


Fic. 2. The variation of the dimensionless wave-length \ 
and lateral range 2max—2min, Of periodic orbits in a linearly 
varying magnetic field, with the maximal extension 2max. 


considered linear with distance to a first approxi- 
mation. The orbits calculated from Eq. (19) will 
therefore represent the actual paths of the ions 
or electrons to at least a first approximation. 
Application of such orbits as in Figs. 1(g) and 
i(f) to magnetic focusing instruments will be 
discussed below. 


EXPONENTIALLY VARYING FIELD 


For the case of an exponentially varying field 
in Cartesian coordinates we may write H as: 


H=Hee*, 
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where Hp is the magnitude of the field strength 
at the value of x chosen as the origin of the 
coordinate system, and b is the constant deter- 
mining the rate of change of H with x. 
The function f then becomes: 
aH» 
where é is a constant of integration which may 


be adjusted to give the orbit the desired slope 
at any value of x. For this case Eq. (8) becomes: 


(29) 


By making the transformation : 


v=——e"*, 
b 


Eq. (29) may be readily integrated to give the 
general solution : 


| (v+2) 


v 
when |2| <1, and 


(30) 


(v+é) 
y v 


v |}: 


when —1. 
For the special case, ¢= —1, 


Orbits defined by Eq. (30) are illustrated by 
(a), (b), and (c) in Fig. 3. In each case the 
particle enters the magnetic field from x= — ~, 
or, physically speaking, from a region of very 
weak H. As must obviously be the case, each 
orbit is symmetrical about the x axis. Each path 
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Y, CENTIMETERS 


(A) (8) 


Nab 


H=H,e* 

H=4 46x10" Oersteds 

V=300 volts 

M= 150 

(c) (0) 


-I0 -8 -6 -4-2 0 


X, CENTIMETERS 


mined by the 
particles of 3 


is shown in the neighborhood of H=H, which 
establishes the origin of the x coordinate. Orbits 
(a), (b), and (c) represent a particle of 150 
atomic mass units and 300 electron volts energy 

.coming from x= — © with initial slope of 1, 0, 
and —1, respectively. Orbit (d) is for a particle 
of the same mass and energy but so projected in 
the x, y plane that the path has infinite slope at 
x=0. 

Orbits of the type illustrated by (b) are of 
particular interest, as the total extension in the 
y direction depends only on the exponential 
coefficient b, and is independent of the values of 
the mass and velocity of the ion or electron. 
This fact is immediately seen upon setting ¢=0 
in Eq. (30), which physically corresponds to the 
particles approaching regions of stronger mag- 
netic field with an original direction parallel to 
the x direction. This situation affords an oppor- 
tunity for perfect focusing, further discussion of 
which will be given later. 

If the ion or electron is introduced into the 
magnetic field in such a manner that its subse- 
quent orbit is bounded in the x direction, it will 
be of the type illustrated by (d), Fig. 3. This 
type of orbit is also of interest as regards applica- 
tion to focusing problems, and will be discussed 
in more detail below. For this type trajectory 
the total width in the x dimension is readily 
calculated from the separation of the roots of 
Eq. (15). 

The behavior of orbits of type (d) as regards 
their total width in the x direction and their 
wave-length \ can be discussed in a manner 


Fic. 3. Various types of orbits possible in an exponential 
rameter ¢, defined by Eq. (28). All orbits are 

-ev energy and 150 atomic mass units. Note: 

explanation at the right of the cut should read Ho=4.46X 10° oersteds. 
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netic field, deter- 


or singly-charged 
The second line of the 


similar to that used for the linearly varying 
field. For a particular value of @ we have: 


: (1—2) 
Xmax 
b 
(32) 


b 
Xm log —-~(—1-2). 
b aH 0 
As ¢<—1 for all such orbits the arguments of 
the logarithms will be positive. We may express 
Xmax—Xmin aS a function of Xmax by eliminating ¢ 
in Eqs. (32) to give: 


2b —), 
Xmax — Xm og 
aH 
which behaves qualitatively like the similar 
quantity for the linear case which is plotted in 
Fig. 2. 

From the nature of Eq. (31) it is seen that the 
orbits are periodic with a wave-length of \ 


given by: 
é 
i+ ). 
b (e?—1)! 


RADIAL FIELDS 


Although Eqs. (10) and (14) make it possible 
to determine the path of a charged particle in 


any 
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any field of circular symmetry, the functional 
yariation of H we shall discuss here is given by: 


where »>—1, and Hp is a constant. It will, of 
course, be physically impossible for H to obey 
this relation for values of r including the origin. 
However, it will suffice if the relation is obeyed 
within certain limits of r, or for values of r above 
a definite radius. 

For this case (and »#1) g(r) becomes: 


aH, é 
+, (34) 
(1—n)r" 


g 


Centumeters 


@ 
(c) +150 
300 electron voits 


Fic. 4. Non-periodic orbits with common asymptotes in 


a radially varying ee field, with n=3, for singly- 
charged particles of ev and various atomic masses M. 


and we have: 


Equation (35) is in general more amenable to 
numerical than analytical treatment. 

One case of important physical interest does, 
however, admit of simple analytical treatment. 
This is that of charged particles approaching 
regions of strong H from regions of weak H, 
and with an initial direction along a radius 
vector and toward the axis of symmetry. 
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Before discussing the solution for this par- 
ticular case it will be of interest to examine the 
general features of the orbits defined by Eq. (35). 
It will be noted first, from Eq. (34), that if 
n>0O (#1), g(r)-0 as Hence by Eq. (13), 
r(d0/dr) +0 as ro. This means that for 
n>0O (#1) all orbits entering the field from 
infinity must be asymptotically parallel to a 
radius vector. As all radii vectors are equivalent, 
these orbits may thus be considered as necessarily 
entering the field parallel to the polar axis 
(@=0). If, in particular, these orbits asymptoti- 
cally degenerate into linear paths at a normal 
distance h from the polar axis, i.e., so that 


rsin 6=h, 
then it follows that in the limit as r>~, 
dé h h 


dr r cos @ r 


By Eq. (34), such an asymptotic behavior can 
be satisfied only if m>1 and ¢= —h. Accordingly, 
in inverse power radial fields for which n>1, —é 
is simply the normal distance from the polar 
axis with which the particle enters the field. 
The restriction 1 evidently means that if 
there can be no orbits which enter the field 
along linear paths parallel to a radius vector. 
If 0<n<1 the value of @ can be determined by 
fixing r(d@/dr) at a given radius, computing g 
from Eq. (13), and then @ by Eq. (34). 

If n=0, the orbits entering or leaving the field 
at r= © will asymptotically cut the radii vectors 
at an angle sin“ aHo. For this case the integral 
for 6(r) can be readily evaluated to give ex- 
pressions similar to Eqs. (30) and (31). If 
—1<n<0 the orbits will be restricted to finite 
regions of the (r, @) plane, and be bounded by a 
maximum value of r. 

Returning to the special case m>1, and also 
assuming that the particle actually enters the 
field along the polar axis (¢=0), we may readily 
integrate Eq. (35) and obtain: 


aH 
+— sin-! ———. (36) 
n (1—n)r” 
Orbits of this type for »=3 are plotted in 
Fig. 4. It is interesting to note that these orbits 
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Fic. 5. Typical orbits in a radially varying magnetic field, 
with »=}3, which are convex to the origin at r=rmin. 


have an angular spread independent of mass and 
energy, quite analogous to those for the expo- 
nentially varying field in Cartesian coordinates, 
shown in Fig. 3(b). These trajectories will be 
discussed in more detail later. 

Other types of orbits obtainable in radial 
fields are illustrated in Figs. 5 and 6. In Fig. 5 
the particles are so projected, in a field with 
n=4, that the paths are convex with respect to 
the origin at the points of closest approach. 
These were determined by numerical integration 
of Eq. (35), after choosing the point of minimum 
approach and determining the corresponding 
value of ¢. In a similar manner were computed 
the orbits shown in Fig. 6, which are concave to 
the origin at the points of closest approach. In 
Fig. 5 may be seen a periodic or cyclic type of 
orbit of the same general properties as possessed 
by the periodic orbits in the Cartesian coordinate 
cases. 

For periodic orbits, the equation 


1—g*=0, 


which is analogous to Eq. (15) for the Cartesian 
case, must have two real and unequal roots. 
For a given value of m and @ it is possible to 
determine whether the resulting orbit will be 
periodic or not by determining whether g will 
pass through both the values —1 and +1 as 7 is 
increased. 
APPLICATIONS 

In addition to illustrating the application of 

the general method of calculation given above to 
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specific cases, the orbits previously discusgeq 
have certain characteristics of promise for instry. 
ments such as mass spectrometers, mass spectro. 
graphs, beta-ray spectrometers, etc. For example, 
Fig. 7 illustrates the orbits of ions of masses 1, 
35, and 150, expressed in atomic mass units, and 
all having a kinetic energy of 300 electron volts, 
The type of equation which describes these 
orbits, as well as that of Fig. 3(b), is obtained by 
setting ¢=0 in Eq. (30), giving 


1 aH, 


From these figures and Eq. (37) it is seen that 
all such ions which initially approach an expo. 
nentially varying field in the direction of ip. 
creasing H will turn and recede in an opposite 
direction. Also, the ions will recede asymptoti- 
cally to a line displaced from their original line 
of approach by an amount independent of the 
momentum of the particle and depending only 
upon the exponential coefficient b. 

With this arrangement all ions beginning at a 
point where the field is negligibly weak, and with 
a direction toward regions of increasing H, will 
return to a single point where the magnetic 
field is also negligibly weak. ; 

A schematic diagram of an instrument which 
might be designed to utilize this type of orbit 
may be seen in Fig. 8. The design shown here is 
for a mass spectrometer, but with obvious 
changes it could be used as a beta-ray spec- 


Fic. 6. Typical orbits in a radially varying magnetic field, 
with n= 4, which are concave to the origin at 
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trometer, or for special studies such as electron 
or ion scattering, etc. By adjusting an aperture 
or slit, located as shown in Fig. 8, the instrument 
may be made to focus either only particles of a 
narrow momentum range, such as ions of a 
definite mass, or particles of a wide momentum 
or energy range. An idea of the amount of 
dispersion taking place at the region of the slit 
may be obtained by inspection of Fig. 7, which 
represents paths of physical interest. The separa- 
tion in the turning points of any two ion beams 
of momentum parameters a; and a2, may be 
readily calculated from Eq. (37), which gives 


a2 
. 


(mv) 


(mv)2 


For an instrument such as this, it would be nec- 
essary for the magnetic field to vary in an expo- 
nential manner only over a region large enough 
to encompass the dimensions of the instrument. 

In designing pole pieces of an electromagnet so 
as to obtain an exponentially varying field, the 
first requirement is to have the dimension in the 
direction perpendicular to the x, y plane long 
compared to the size of the instrument desired, 
so as to allow the field to be accurately two 
dimensional. Next it would be necessary that 
the design be such that the three equations: 


=(@, (38) 
0x? 92? 
Ox 
-at z=0, (39) 
—= Hye 
Oz J 


be satisfied, where M represents magnetostatic 
potential. 
A solution of Eq. (38), which also satisfies 
Eqs. (39), is: 
Hye™* 


sin bz. 


M= 


If the magnetostatic potentials of the opposing 
pole faces are +My and — Mg, respectively, we 
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Fic. 7. Typical orbits with common asymptotes in an 
exponentially varying magnetic field. 


have the equation : 


1 bM 
b Ho 


(40) 


to define the contours of the pole faces. 

In Fig. 4 are seen orbits which, for the case of 
radially symmetric fields, are analogous to the 
orbits just discussed. They are defined by Eq. 
(36). It is obvious that since these orbits are 
analogous to those of Figs. 7 and 5, the same 
focusing properties could be applied in corre- 
sponding instruments. As for the Cartesian expo- 
nential field, the optics of the radial field system 
is such that ions originating at a point where 
the field is negligibly weak, i.e., at “infinity,” and 
with an initial velocity toward the center of 
symmetry of the field, will be all focused, regard- 
less of their individual momenta, at another 
single point. This point will be on a radius 
vector which is rotated from the line of approach 
radius vector by an angle which depends only 
upon the exponent in the equation defining the 
variation of H with r. It must also be remembered 
that orbits with this property are possible only 
for n>1 in Eq. (34). 
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_ A schematic diagram of an instrument de- 
signed to employ this type of orbit is shown in 
Fig. 9. An idea of the dispersion obtainable may 
be seen from Fig. 4, and may be expressed 
analytically as: 


Ar 1Aa 1 A(mv) 


where r and r+Ar are the turning radii for 
particles of parameters a and a+Aa. Fields which 
decrease with radial distance with the same 
functional form have been successfully used in 
the betatron.? 
To obtain the field variation just discussed the 
| circularly symmetric pole faces must be so 
shaped that the following equations are satisfied : 


aM aM 
—(sina—)=0, (41) 
or sin & 


(42) 
-(— Ho / 


| ? D. W. Kerst, Phys. Rev. 60, 47 (1941). 


of a magnetic analyzer utilizing an 
exponentially varying magnetic field. 


_ from one pole face to the other, so as to ensure 


where M is the magnetostatic potential, and 4 jg 
the latitude angle. 

Setting and M=(H)/r*) F(x), we find 
F(u) must satisfy the Legendre equation : 


d?F adF 
(u? — 1) — n(n—1) F=0. 
du 


To fix the specific solutions for F(u) applying 
to the present problem, the equivalents of Eg. 
(42), namely: F(0)=0; F’(0)=—1, are to be 
imposed on the functions F(y). 

If we let the opposite pole faces have magneto- 
static potential +Mp> and — Mo, the contours 
will be given by: 


(43) 


the sign of F(@) automatically changing in passing 


a real value for r. 

As the region over which it is desired that 7 
have the specified functional form will not extend 
to the origin, the pole faces in this region need 
not follow the shape required by Eq. (43). This 
physical observation also obviates the apparent 
difficulty arising from the singularities at .=+1 
in the functions F(x) when n is not an integer. 


Fic. 9. 
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The periodic or cyclic orbits of the type in 
Figs. 1(f) and (g), and Fig. 3(d), also have 
properties of promise for focusing applications. 
One of the most important of these properties 
js illustrated in Fig. 2. This shows how the 
wave-length X, or spatial periodicity of the orbit, 
decreases with increasing Xmax- This curve was 
drawn for H varying linearly with x, but its 
character will be essentially the same for other 
functional variations wherein H increases mono- 
tonically with increasing x. 

An application of such orbits may be made to 
the design of a mass spectrometer or similar 
instrument, as indicated schematically in Fig. 10. 
In such an instrument the ions are so acceler- 
ated by a suitable source that they are projected 
into the magnetic field with the xmax for each 


orbit lying within the exit slits. 


In Fig. 11 may be seen a schematic diagram of 
the state of affairs at the exit slits S,; and S, of 
such an ion source. We shall assume that the 
potential distribution in the source is such that 
all ions acquire essentially the same potential, 
and that the last two slits, ie., S; and S, in 
Fig. 11, are at the same potential and are 
operative only in defining the emergent beam. 
Moreover, these slits will be assumed to be 
geometrically similar and located identically as 
regards the x coordinate. 


Fic. 9. Schematic diagram of a magnetic analyzer utilizing 
a radially varying magnetic field. 
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Gas Inlet 


To Vacuum 


Fic. 10. Schematic diagram of a magnetic analyzer 
utilizing _—— orbits in an exponentially varying mag- 
netic field. 

With such an arrangement all ions leaving the 
exit slits will follow orbits having approximately 
the same radius of curvature at the slits, differing 
only in the angle of emergence. 

‘In Fig. 11, (c) represents the orbit with the 
smallest value of xmax that can emerge from the 
slits, (d) that with the largest value of xmax, 
while (a) and (b) represent the paths of greatest 
divergence in angle of emergence. As S; and S; 
cover the same range along the x direction, the 
Xmax Will be the same for orbits (a) and (b). 

Consider now the orbits (a), (b), (c), and (d) 
when the particles have advanced by exactly 
one wave-length in their periodic paths. Since 
(a) and (b) have the same Xmax, they will have 
the same \. Therefore a slit S;’ of the same width 
as S; and located at the same values of x, but 
displaced along y by a distance X, will collect 
both (a) and (b). That is, (a) and (b) will pass 
through S,’ with exactly the same geometry as 
they passed through S;. As orbit (c) has a smaller 
value of xmax, it will be displaced slightly upward, 
relative to (a) and (b). Conversely, as orbit (d) 
has a larger value of xmax, its \ will be smaller 
than that of (a) and (b). Hence as it passes 
through S;,’ it will be displaced slightly downward 
relative to (a) and (b). 

This shifting of the orbits (c) and (d) relative 
to (a) and (b) means that the latter will limit the 
lateral spread in the ion beam when entering 
S,’ just as they defined the angular spread in 
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Fic. 11. Schematic dia of the geometry of the 
orbits at the source and collector slits of an instrument of 


the type shown in Fig. 10. 


the beam emerging from the source slits S; 
and S:. In the design of such an instrument, as 
schematically shown in Fig. 10, the final two 
exit slits will correspond to S; and S2, and the 
collector slit should be of the same dimensions 
as S; and located at the same distance along the 
x axis. The separation along the y axis between 
S, and S;’ will then define the wave-length \ for 
which the above-described focusing process 
applies. Any ion may be made to follow an 
orbit of any chosen \ by properly fixing the 
accelerating voltage in the ion source. It is to be 
noted that for the collector the one slit S,’ is 
sufficient. S2’ is shown in Fig. 11 only to indicate 
the geometry of the orbits at this point. If the 
instrument is ‘to be used so that only ions of 
one mass or of small range of momentum are to 
pass the collector, an aperture or slit may be 
introduced, as shown in Fig. 10. 


The calculation of the dispersion at the region 
of the selector slit or aperture is easily made by 
application of Eq. (29). For this purpose it js 
convenient to choose the origin at the exit slits 
of the ion source. Then, since dy/dx= © at the 
ion source (x =0), @ must have the value 


aH 0 


At the turning point of the orbit, i.e., at the 
plane of the selecting slits, dy/dx= © again, so 
that: 


—1, 


giving for the location of the turning point: 
x=- —-—+}. 
b aH 


For two ion beams of momentum parameters a, 
and a2 the turning points will then be separated 
by a distance: 


1 


A similar type of instrument could be designed 
using the periodic orbits in a radial field, such 
as shown in Fig. 5. Similar considerations with 
respect to focusing and the use of selector slits 
will be applicable to such instruments. Finally, 
it should be noted that the above discussion has 
referred to particles moving in the median plane 
between the pole pieces generating the magnetic 
field. Hence in practical designs account will 
have to be taken of the perturbations in the 
orbits of particles off the median plane owing 
to the curvature of the magnetic flux lines. 

The authors are indebted to Mr. N. F. Kerr 
for assistance in preparing the drawings and 
making the calculations, and to Dr. P. D. Foote, 
Executive Vice President of The Gulf Research 
and Development Company, for permission to 
publish this paper. 
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The X-Ray K-Absorption Edge for Sodium and the Secondary Structure* 


Accompanying the Edge for the Sodium Halides 


KENNETH C. RULEt 
Institute of Physics, University of Upsala, Upsala, Sweden 
(Received May 24, 1944) 


A high vacuum, curved-crystal spectrograph was constructed for study of the K-absorption 
edge and the accompanying secondary absorption structures for evaporated films of sodium 
metal and of the sodium halides. Comparison of the K8, emission line and the K-absorption 
edge for metallic sodium indicates that the K8, line is not semi-optical in origin, but represents 
a transition from the My 1m valence shell. The energies of the 1s, 2p, and 3p electron energy 
levels for metallic sodium are calculated. The energy difference between the Kf; emission line 
and the K-absorption edge for Na* in NaCl is found to be 5 ev greater than the difference 
between the K@; emission line and the K-absorption edge for Cl- in NaCl. Strong secondary 
absorption patterns are found on the short wave-length side of the K-absorption edge for Na* 
in the sodium halides. The distances from the various maxima and minima in the pattern to 
the main edges are correlated with the lattice distances and satisfactory agreement with 
Kronig’s hypothesis is found except for the features near the main edge in the spectrum for 
NaF. Fair agreement is also found in comparison of these patterns with the patterns accom- 
panying the K-absorption edge for K* in the potassium halides, as determined by Brewington. 


N Kronig’s' hypothesis of the origin of the 

secondary absorption structure which often 
appears on the short wave-length side of x-ray 
absorption edges of elements in crystal lattices, 
it is supposed that the final energies of the elec- 
trons ejected in any given direction through the 
lattice are restricted by the occurrence of Bragg- 
like reflections in the periodic potential field in 
the lattice. Owing to the overlapping of the effects 
of reflections from numerous sets of planes, actual 
calculations of secondary absorption patterns 
have not been carried out rigorously even for the 
simplest lattice types, although approximate 
calculations have been attempted.2 A number of 
simple deductions from the theory have, however, 
been subjected to extensive experimental tests in 
which metals and binary alloys have been em- 
ployed as. absorbing media. Most of these ex- 
perimental data have been found to be harmony 
with the theory,’ but in the case of alloys the 


* This work was carried out at the Institute of Physics 
of the University of Upsala in Sweden and was aided by 
the G. Hilmer Lundbeck fellowship of the American- 
Scandinavian Foundation. 

tAuthor’s present address: South Charleston Division, 
Westvaco Chlorine Products Corporation, South Charles- 
ton, West Virginia. 

'R. de L. ee Zeits. f. Physik 70, 317 (1931); 
rm on and J. Veldkamp, Zeits. f. Physik 74, 191 

*R. de L. Kronig, Zeits. f. Physik 75, 191 (1932); 
C. — J. de phys. et rad. 1, 133 (1940). 

*D. Coster and H. Levy, Physica 6, 44 (1939). 


crystal lattice is not always the decisive factor in 
determining the secondary structure. 

The most serious discrepancies with theory 
have been encountered in the case of absorption 
edges of elements in ionic lattices.*-* Since 
Kronig’s theory applies only to the motion of a 
“free electron,’ discrepancies in the structure 
located within twenty-five or fifty electron volts 
from the edge can be attributed to effects of the 
field about the ion from whose x-ray level the 
electron is ejected. Stephenson® has found, how- 
ever, that the ionic charge and the ionic group 
from which the electron is ejected affect the 
secondary structure at least as far as 100 ev from 
the edge. 

It has been deduced from Kronig’s theory that 
for an ion occurring in a succession of ionic 
lattices of a given type and varying lattice dis- 
tances, the difference in energy between a given 
feature in the secondary absorption structure and 
the main edge should be inversely proportional to 
the square of the lattice distance. Since the 


‘J. H. Munier, J. A. Bearden, and C. H. Shaw, Phys. 
Rev. 58, 537 (1940). 
asin” Bearden and H. Friedman, Phys. Rev. 58, 387 
*Vola P. Barton and Geo. A. Lindsay, hy Rev. 46, . 
362 (1934); S. Yoshida, Sci. Papers Tokyo Inst. Phys.- 
Chem. Research 38, 272 (1941). 
7 D. Coster and G. H. Klamer, Physica 1, 145 (1933). 
8G. P. Brewington, Phys. Rev. 46, 861 (1934). 
*S. I. Stephenson, Phys. Rev. 58, 873 (1941). 
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absolute separation depends only on the lattice 
and is independent of the energy of the main 
edge, it is evident that the relative separation of 
the secondary absorption structure from the 
main edge should be greatest’ for an edge of long 


SSS 


MAAN 


Sem 
Fic. 1. Cross section of the cathode-ray oven. 


wave-length. Previous experiments with ionic 
lattices as absorbing media have been confined to 
K and L edges of comparatively short wave- 
lengths. Early indications were that secondary 
absorption becomes weaker for elements of low 
atomic number,’ but recent work has succeeded 
in bringing out the secondary structure in the 
K edges for metallic Al and Mg.‘ In view of the 
realizable gain in resolving power, it therefore 
seemed worth while to attempt to study the 
structure accompanying an edge at still longer 
wave-length, such as the K edge for Na in its 
halides. A second advantage in a study of the 
sodium halides is that the lattice constants in- 
crease as much as 1.4-fold from the fluoride to the 
iodide, a range which should make it possible to 
draw significant conclusions as to the correctness 
of Kronig’s hypothesis in respect to dependence 


_ of the separation from the main edge on lattice 


distance. 


Aside from the objective of studying the 


secondary absorption structure, a determination 


of the wave-length of the main sodium edge, ip 
which the various halides and also sodium metal 
were used as absorption films, was of interest jn 
itself. The only previous determination found in 
the literature was an incidental measurement by 
Johnson,'® who employed sodium chloride as the 
absorption film. 


EXPERIMENTAL 
Apparatus 


The apparatus and technique employed were 
designed to avoid excessively long times of ex- 
posure in bringing out absorption edges and 
secondary absorption structure in the region of 
long wave-lengths. The curved-crystal type of 
spectrograph was therefore chosen. A further 
advantage of this instrument is that the edge and 
the absorption spectrum extending far out from 
the edge can be photographed on a single plate at 
one setting of the crystal. The high vacuum ex. 
tended throughout the spectrograph so that no 
diaphragm was required over the mouth of the 
x-ray tube. Since thin, homogeneous absorption 
films would further contribute to lowering the 
exposure time, the films were formed by con- 
densation on an aluminum foil of vapors produced 
by heating the materials in a cathode-ray oven. 
In order that films of metallic sodium as well as 
of the sodium halides could be deposited and 
employed without exposure to air, the cathode- 
ray oven was constructed as an integral part of 
the high vacuum spectrograph. 

The cathode-ray oven, Fig. 1, was essentially a 
water-cooled, metal, x-ray tube mounted at the 
bottom of the crystal chamber. The anticathode 
was designed" to minimize loss of heat by radia- 
tion and was thimble-shaped in order to facilitate 
heat transfer to the material to be evaporated. A 
loop of tungsten wire at ground potential was the 
source of an electron beam of 1 to 200 milliamperes 
which served to maintain the crucible at 600 to 
1400°C when a potential of 1 to 20 kv was applied 
to the anticathode. The oven was evacuated 
through a flexible metal tube connected directly 
to the oil-diffusion pump. ae 

The x-ray tube was a water-cooled, metal tube 
which projected through the cover of the crystal 


” Nils G. Johnson, Dissertation (Lund, 1939), p. 50. 
"H. M. O'Bryan, Rev. Sci. Inst. 5, 125 (1934). 
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chamber, Fig. 2. Since a bent crystal requires a 
long focal spot, the cathode was a straight section 
of tungsten wire screened so as to produce a focal 
spot 18 mm long and 1 mm high. The aluminum 
foil, upon which the vapors were condensed to 
form the absorption film, was mounted so that it 
could be revolved about a vertical axis and stand 
either beside the mouth of the oven or in the path 
of the beam from the x-ray tube to the crystal. 
When the film was beside the mouth of the oven, 
it could be lowered over the oven by means of an 
arm attached to a ground metal joint mounted 
on the side of the crystal chamber. All other 
joints except the two for the manipulation of the 
film were flanged joints sealed by rubber-ring 
gaskets mounted in a groove in one of each pair of 
flanges. 

The crystal holder and the cover carrying the 
x-ray tube could be rotated independently to give 
any glancing angle from 15° to 80°. The holder 
for the bent crystal was of the type designed by 
Haglund and described by Ohlin,” the opening in 
the front face being 1.5 mm by 20 mm. 

The beam reflected by the crystal passed 
through a two-inch hole drilled through the 
cylindrical wall of the crystal chamber and 
through a cast-iron plate two inches thick and 
one foot square, to which the wall of the chamber 
was fastened. A shelf of cantilever construction 
was bolted to the other side of the cast-iron plate 
in order to support a track upon which the table 
bearing the plate holder moved. In order to 
protect the plate from visible light and scattered 
radiation, the plate holder was surrounded by an 
aluminum shield in which a window covered by 
two thicknesses of 0.64 aluminum foil was pro- 
vided for passage of the x-rays. 

A flanged, steel pipe three feet long, ten inches 
in diameter, and blinded off at one end was 
mounted on a track to form the vacuum chamber 
for the plate holder when rolled forward to abut 
against the cast-iron plate. The degree of evacua- 
tion at any time could be determined by means 
of a Pirani tube and an ionization gauge, both 
permanently attached to the crystal chamber. 
The measuring circuit for these gauges was built 
according to the design of Sten v. Friesen." 

® Per Ohlin, Dissertation (Upsala, 1941), pp. 25, 26. 


*Sten v. Friesen, Arkiv. f. Mat. Astronom. Fysik B27, 
No. 11 (1939). 


Operation 

The process of forming the absorption film was 
not started until the spectrograph and cathode- 
ray oven had been evacuated to 2X10~ mm of 
Hg or less. To produce films of sodium metal a 
reaction mixture of equal parts by weight of 
powdered zirconium and sodium molybdate was 
employed.“ The films of sodium halides were 
produced by evaporating C. P. halides. The 
evaporation process was always carried out by 
increasing the strength of the cathode-ray current 
in small steps, with sufficient time between suc- 
cessive steps to reduce each temporary rise in 
pressure. Whether or not a deposit had been 
evenly formed on the aluminum foil could be 
determined by tilting up the film holder until the 
under side of the foil could be viewed through a 
glass window in the side wall of the crystal cham- 
ber. The approximate thickness of the deposit 
was determined, after recording the absorption 
spectrum, by weighing the film and the aluminum 
foil. 

A tungsten anticathode was employed as the 
source of the continuous radiation since it pro- 
duced ‘no interfering lines in the region of the 
sodium edge. Because the rectifying equipment 
available, four tubes in Greinacher coupling 
without condensers, did not prevent a half-sine- 
wave fluctuation of voltage, the effective voltage 
had to be maintained at or below 1.6 kv to pre- 
vent higher orders of radiation from obscuring 
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Fic. 2. Cross section of the crystal chamber. 


™“ The use of this ratio was an oversight, since de Boer, 
Broos, and Emmens [Zeits. f. anorg. allgem. Chemie 191, 
113 (1930) ] state that the product evaporated from a 1 to 1 
mixture is sodium metal containing some oxide, while a 
mixture of 1 part sodium molybdate to 4 parts zirconium 
gives the pure metal. 
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the absorption spectrum at 1.1 kv. The Ka; line 
of sulfur always appeared on the plates, but not 
the Ka; line of potassium, which indicated that 
some second-order, but no third-order radiation 
was excited. Agfa Isochrome Plates, 18/10 Din, 
were employed for recording the spectra although 
their large grain size set a low limit to the 
enlargement which could be employed in making 
photometer curves. 

Reference lines were exposed on one-half of the 
photographic plate either before or after the 
absorption spectrum was recorded. KCl and ZnO 
powders on a copper base on the other half of the 
anticathode were the source of the reference lines, 
19-K Ka, at 3733.68 x.u. (effectively 11217.05 x.u. 
for a mica crystal in the third order) and 
30-Zn Kay, at 1432.17 x.u. (11474.56 x.u. for mica 
in the eighth order).'* 


Measurements and Errors 


A scratch was made on each plate parallel to 
the absorption edge, and the distances from this 
scratch to the reference lines were measured on a 
comparator. All measurements of positions of the 
edges and secondary absorption structures were 
made on photometer curves with respect to the 
sharp peak produced by this scratch. An abscissa 
enlargement of three was found suitable, the 
exact enlargement always being determined by 
recording at the same time a photometer curve of 
two scratches on a plate, the distances between 
the scratches having been determined by means 
of the comparator. A millimeter netting was also 
copied onto the photometer curve in order to 
provide automatic correction for shrinkage of the 
paper during photographic development. 

The position of the edge and other main fea- 
tures on the photometer curves could be measured 
to +0.1 mm, corresponding to +0.03 mm on the 


plate or to +1.1 x.u. The positions of the refer- 


ence lines could, of course, be determined with 
greater precision, while the precision was much 
lower for the weaker features of the spectra. 

A number of systematic errors in wave-length 
measurements arise when a bent-crystal spectro- 
graph is employed and the spectrum is recorded 


on a plane plate tangent to the focusing circle. 


% Manne bahn, Spektroscopie der Rontgenstrahlen 
(Julius ten Berlin, 1931). 
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Equations for calculating the various errors have 
been published by Sandstrom.'® These errors 
arise from the geometry of the system evep 
though the curved surface of the crystal is per. 
fectly cylindrical. For the spectrograph employed 
in this work, the geometrical data required for 
Sandstrom’s equations are: radius of the bent 
crystal 498 mm, crystal opening 20 mm X1.5 mm, 
height of focal spot 1 mm, glancing angle 35° 49’, 
distance from focal spot to crystal 80 mm, and 
thickness of the photographic emulsion 12y, 
Variations in temperature of the crystal are also 
a factor, and while the temperature of the crystal 
itself was not measured, the maximum tempera- 

ture variation in the room during an exposure 
was 3°C. 

At a crystal opening of 20 mm, which was that 
employed in obtaining most of the data, the 
largest error under these conditions is the 
focusing error, which arises because of the fact 
that the radius of curvature of the surface of the 
bent crystal is twice that of the focusing circle. 
When a crystal opening of 20 mm is employed, 
the focusing defect at a glancing angle of 35° 
should broaden a sharp line toward shorter wave- 
lengths into a beam 4.7 x.u. wide. A reduction of 
the crystal opening to 8 mm would reduce the 
focusing defect to 1.1 x.u., which should be 
practically counterbalanced by other effects so as 
to produce symmetrical broadening. Comparison 
of the positions of the reference lines when 
recorded with 20-mm and 8-mm openings re- 
vealed a shift toward shorter wave-lengths of 
only 0.4 x.u. at the larger opening. 

In addition to errors owing to geometry of the 
system, the possibility of lack of definition owing 
to imperfect bending of the crystal and inaccurate 
focusing of the apparatus must be considered. 
The radius of curvature of the crystal was de- 
termined by the standard methods employed at 
the Institute of Physics at Upsala, which involve 
observation of the points of focus of x-ray emis- 
sion lines on a nearly horizontally inclined photo- 
graphic plate. From the radius of curvature so 
determined the required distance from the crystal 
to the plate can be calculated for any glancing 
angle. In order to check the sharpness of focus, 


16 Arne Sandstrom, Zeits. f. Sag 84, 541 (iss 
622 (1934); Dissertation (Upsa ” 1935), BP. 21-27; also 
H. H. Johann, Zeits. f. Phys 69, 185 (1931 
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TABLE I. Absorption films, exposure times, and individual measurements of the main edges. 


Wave-lengths for the main edges 
x.u. 


Ki Ke 


S885 


NaBr 0.3 53 15 20 53 
0.5 54 8 20 54 
20 55 


Nal 


0.6 56 16 
0.2 59 22 


59 


11544.3 
11544.6 
11545.5 
11545.5 


11488.8 
11480.6 
11484.5 


11509.4 
11507.1 
11506.5 


11523.8 
11520.1 
11522.1 


11498.9 
11502.3 
11498.9 
11501.2 
11500.6 
11503.1 
11503.1 


11500.3 
11502.6 
11501.2 
11502.3 


11504.6 
11505.1 
11504.6 


11507.1 
11502.9 


11517.5 
11518.5 
11519.8 
11519.0 
11522.0 
11523.7 
\ 


* Weight given for oxidized film after removal from spectrograph. 


estimates were made of the widths of emission 
lines exposed at large glancing angles with the full 
opening of the crystal. At a glancing angle of 58° 
the full width for 29-Cu Ka; (eleventh order) was 
found to be 0.6 x.u. at half-maximum intensity, 
which is comparable with the widths observed in 
a double-crystal spectrometer.!” It is evident that 
even with the full crystal opening no significant 
broadening occurred such as might be expected 
to result from imperfect bending of the crystal. 
It can be concluded that the systematic errors 


TaBLe II. Average wave-length and energy values for 
the sodium K-absorption edge. 


Wave-length 
Absorbing Average Energy values 

material Edge x.u. deviation v/R ev Aev 

Na Ki 11545. 0.6 78.932 1071.48 
NaF K’ 11520. 79.103 1073.79 2.32 
Ki 11501.2 lis 79.233 1075.56 4.09 
NaCl Ki 11501. O.» 79.225 1075.52 4.05 
NaBr 1 11504.s 0.3 79.208 1075.22 3.75 
Nal Ki 1150S.0 2a 79.206 1075.20 3.73 
Ki 11522. lia 79.090 1073.62 2.15 
2 11507.6 0.9 79.189 1074.96 3.49 
K; 11484.6 2.8 79.347 1077.11 5.64 


Lyman G. Parratt, Phys. Rev. 50, 1-15 (1936). 


are smaller than the accidental errors in de- 
termining the wave-lengths. However, repetition 
of the measurements at higher dispersion might 
reveal additional structure in the edges, in which 
case some of the wave-lengths obtained might 
differ appreciably from those reported here. 


DATA 


Table I lists each absorption film from which 
absorption spectra were obtained, together with 
the wave-length values found for each of the 
main edges. Table II gives the average value 


_ obtained for the wave-length of each edge and the 


corresponding average energy values, calcu- 
lated from the relations v/R=10"/109737\ and 
ev=12370.2X10°/A, where A is in x units. 
Average wave-lengths for the edges and all re- 
producible features in the secondary absorption 
structures are listed in Table III. The figures in 
the third column of each compound in Table III 
indicate the average deviation of the individual 
measurements from the average, and the figures 
in the fourth column indicate the number of 
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TaBLeE III. Wave-length values for the features of 
the sodium K-absorption edge. 


Z 
2, 


No. of 
Av. Av. measure- Av. A 
Feature x.u. dev. ments Feature x.u. 


NaCl 
kK’ 
Ki 


oo 


in 


K 
A 
a 
B 
b 
8 
D 


* Film of metallic sodium oxidized by exposure to air. 


measurements (number of plates) entering into 
the determination of each average. Typical 


Fic. 3. Photometer curves of the for a film of sodium 
metal and for an air-oxidized film. 


Fic. 4. Photometer curves of the edges for the sodium 
halides. Reading from top to bottom, the curves corre- 
spond to plates Nos. 59, 54, 28, and 49. 


photometer curves are reproduced in Figs. 3 
and 4. 

In all tables the main edge is distinguished by 
K,, the minima on the photometer curves (i-., 
maxima of absorption) are symbolized by Roman 
capitals, and the maxima on the photometer 
curves (minima of-absorption) by Greek letters. 


\| 


wn 


fi 
‘ 
if 
1481.2 2.9 (7) 11479.0 > Vir 
114264 1.0 (7) 11424.6 4) 
ii 113991 (7 11412.1 4) Ae 
113712 09 (3 11402.1 4) we 
113531 16 (3) 11383.5 1.2 (4) 8 
113475 5.5 (4) 113236 3.6 (4) F 
113327 3.2 (4) 112863 4.3 (4) 
11299.7 43 (7) 112169 60 
112236 3.2 111645 9.0 (4 
111569 70 (7 1100945 5.2 (4) 
10869.1 8.7 (7) | MA 
1 115048 (3) (3) 
1149554 OS (3 2) 
114796 19 (2 2) ‘A 
114218 1.1 (3) (2) 
114128 0.9 (3) (3) 
113995 14 (3 C —-:11093.7 (1) 
| 113358 53 (3) 10932.8 (1) 
11298.2 7.9 (3) 
112448 33 (2) 
111985 9.7 (3) 
11032.7 (1) 
109458 17 (2) | 
Nal Na : - G 
115050 (2.4 (2) Ki 115450 06 (4) 
"414990 (2) K: 115164 74 (4) 
114878 0.6 (2) Ks 114894 08 (4) 
i, 11430.1 B 10958.2 1 
11420.8 (1) C  10677.6 (1) Poe 
co 
t 
oF 
A 
| 


TABLE IV. The energy levels for sodium metal. * 


Level ev 
K 1s 1071.5 78.932 
Lum 2p 30.7 2.269 
Mui 3p 0.64 0.045 


=a 


These symbols have been written in on the 

hotometer curves. The small Roman letters in 
Table II1 designate less conspicuous features, 
whose positions are not indicated in Figs. 3 and 4, 
but which can perhaps be discerned in relation to 
the designated maxima and minima. 

In connection with the discussion to follow it 
should be noted that only those features which 
are similar in appearance (width, strength, and 
shape) as well asin relative position have received 
like designations in the photometer curves of the 
spectra for the several halides. The secondary 
structures for the spectra of sodium metal and 
sodium oxide are quite distinct from those of the 
halides, however. The features a, B, 8, and C in 
the spectrum for sodium metal and C and 7 in 
the spectrum for sodium oxide appeared in each 
case on only one plate and may be fictitious. 


DISCUSSION 
(a) The Main Edges 


The K8, emission line for Na and Mg has been 
considered in the past'*® as representing a semi- 
optical transition since this line is assigned to the 
transition K—My im, and there should be no 
electrons in the My m shell for these elements. 
Wetterblad’® plotted the wave-length differences 
between the Kay. and the K£, emission lines for 
the elements from-17-Cl to 11-Na and found that 
the curve changed sharply in slope at 12-Mg and 
11-Na. A plot of the difference between the K 
edge and the K@, line in terms of ev, Fig. 5, gives a 
smooth curve from Si through Na, however. The 
energies of the K-absorption edges also agree 
within 1 or 2 ev with the energies of the K@, lines. 
Since the K8, emission actually constitutes a 
band, careful determination of the short wave- 
length edges of these bands for the elements from 
17-Cl to 11-Na would undoubtedly result in a 


* E. Biicklin, M. Siegbahn, R. Thoraeus, Phil. Mag. 49, 
513 and 1326 (1925). Referenee 15, p. 392. 
* Thorsten Wetterblad, Zeits. f. Physik 42, 603 (1927). 


K-ABSORPTION EDGE FOR SODIUM 


205 


closer agreement with the K-absorption edges. It 
therefore appears unnecessary to consider the 
Kf, line as a semi-optical transition for Na and 
Mg, but as a transition from a valence shell, 
which may be designated My m. Farineau*® has 
previously reached a similar conclusion as to the 
origin of this line in the case of Mg, Al, Si. With 
this assignment for the K8, emission line, the 
wave-lengths for the K-absorption edge and the 
KB, emission line," together with the value of 
11.9093A = 11885.2 x.u. obtained by Tyrén® for 
the Ka; line, can be used to calculate the energies 
of the 1s, 2p, and 3p levels for metallic sodium. 
The energy values obtained are shown in 
Table IV. 

The main edge for sodium metal is followed by 
two very weak edges, K: and K3. Their wave- 
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Fic. 5. Energy differences for K8,— Ka; and K —Kp,, 
for the elements sodium 
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lengths correspond approximately to the average 


_ of the wave-lengths of the K; and K» edges and 


to the wave-length of the K; edge of the oxide, 
respectively. As mentioned in the experimental 
section, there is reason to believe that the sodium 
film contained some oxide. 

The main edge for the oxide lies 23 x.u. toward 
shorter wave-length, or 2.1; ev toward higher 
energies, than the edge for the metal. Further- 
more, it no longer coincides with the K8, emission 
line for the oxide, which is shifted 43 x.u. in the 
opposite direction as compared with the K@, line 
for the metal." The directions of these shifts are 
the same as are found in analogous cases for 
elements in higher series of the periodic system.” 

The K-absorption edge for Nat in its halides 
is shifted still farther than for the oxide toward 
shorter wave-lengths. The shift is some 40 x.u., 
or about 4 ev toward higher energies. From NaF 
to Nal the shift decreases from 4.09 to 3.73 ev, a 
shift which exceeds the probable error of the 
data. In terms of energy the shifts are somewhat 
larger than the corresponding shifts observed by 
Lindh™ for the K-absorption edge for K* in its 


23V. Hugo Sanner, Dissertation (Upsala, 1941), p. 93. 


* Axel E. Lindh, Arkiv. f. Mat. Astronom. Fysik 18, 
No. 14 (1924). See also reference 15, pp. 289-290. 


halides. Our result for the K edge for Na* jn 
NaCl can be compared with the value obtained 
by Johnson”® in the course of a study of the effects 


- of chemical combination on emission spectra. He 


found a wave-length of 11498 x.u., compared 
with the present result of 11501.6 x.u. In terms of 
energy the difference is 0.4 ev. Since Johnson's 
measurement was not intended to be a precision 
determination, the agreement between the two 
results is satisfactory. 

In Table V the wave-lengths of the absorp- 
tion edges for the halides are compared with 
Wetterblad’s measurements of the KA; line for 
Nat in its halides. The energy for the K8; emis- 
sion line for Nat+ in NaCl is 12.1 ev less than the 
energy for the absorption edge. Valasek®* has 
plotted the separations in units of ev between the 
various K-emission lines and the K-absorption 
edge for Cl- in NaCl. He predicted that if the 
difference between the K@, line and the K-absorp- 
tion edge for Nat in NaCl were plotted in the 
same figure, the KA; line for Nat would coincide 
with the K@; line for Cl-, which lies 7.1 ev from 
the edge. The difference of 5.0 ev between the 
predicted and observed positions for the K8; line 
for Nat is fully an order of magnitude greater 
than the probable experimental errors. It seems 
necessary to conclude either that the K@, 
line for Na+ in NaCl does not represent the 
Nat+—is—Cl-—3p transition as predicted, or 
that the end levels for the K-absorption transi- 
tions for Na+ and Cl- do not both begin at the 
Cl- continuum as assumed. Different end levels 
might be attributed to different selection rules.’ 


% Joseph Valasek, Phys. Rev. 53, 274 (1938). 


TaBLeE VI. Distance of the secondary-structure features to the main edge, in x.u., ev, and a* (ev). 
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Fic. 6. Separation in electron volts between the secondary 


structure and 


the main edges for Na* and K*. + indicates data from Coster and Klamer 


In the case of KCI, also, the coincidence between 
the K@; line for K+ and the K@, line for Cl- is not 
perfect, the K@;5 line for K* lying 1.2 ev farther 
out from the edge. Valasek attributed this dis- 
crepancy, however, to an error in Stelling’s 
energy value” for the K-absorption edge for Cl- 
in KCl. 


(b) Secondary Absorption Structure 
for the Halides 


As was mentioned in the introductory section, 
in view of the difficulties of carrying out detailed 
calculations of the shape of the secondary ab- 
sorption pattern which should arise from the 
overlapping of the numerous forbidden energy 
zones predicted by Kronig’s theory, the theory 
has usually been indirectly tested in case of ionic 
lattices by comparison of the secondary struc- 
tures for a particular absorption edge for a series 
of similar lattices of varying lattice distances. 
For a cubic type of lattice the separation of the 
structure from the edge in terms of energy should 
be inversely proportional to the square of the 
lattice distance. Brewington® applied this crite- 
rion in the case of the secondary absorption 
patterns accompanying the K-absorption edges 
for K+ and Cl- in the series of potassium halides 


* Reference 15, p. 282. 


(reference 7). Other data for K halides taken from Brewington (reference 8). 


and alkali chlorides. He came to the conclusion 
that his data did not agree with the hypothesis. 
The same criterion will now be applied to the 
data given in Table III, and since the potassium 
and sodium halides belong to the same crystal 
type and the K-absorption edges for sodium and 
potassium should represent the same transition, 
the data in Table III will be compared directly 
with the available data for the K edge of K* in 
the potassium halides. 

The separation of each feature of the secondary 
structure from the main edge has been calculated 
from Table III and is listed in Table VI for the 
absorption spectrum from each halide. The sepa- 
rations are expressed in x.u., electron volts, and 
in electron volts multiplied by the square of the 
lattice constant of the halide. The values em- 
ployed for the squares of the lattice constants are 
21.34, 31.65, 35.52, and 41.73 cm X10~"* for NaF, 
NaCl, NaBr, and Nal, respectively. In Fig. 6, the 
separations shown in Table VI in terms of elec- 
tron volts are plotted together with similar data 
of Brewington® and of Coster and Klamer’ for the 
secondary absorption structure accompanying 
the K-absorption edge for potassium in its 
halides. This comparison does not reveal very 
striking likenesses between the spectra for the 
different halides, except for the features close to 
the main edges. In Fig. 7, however, in which the 
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abscissas represent separations multiplied by the 
squares of the lattice constants, the similarities 
are more marked. It is true that the positions of 
the first two maxima and minima for the other 
halides are not matched in the spectra for NaF 
and KF. On the other hand, considering the 
spectra for the sodium halides alone, the double 
minimum E-F is exceptionally strong for NaF, 
obvious for NaCl, and still apparent, though very 
weak, for NaBr, and they lie on a nearly vertical 
line in Fig. 7. Aside from the discrepancies men- 
tioned for NaF and KF, even the differences be- 
tween the patterns for the sodium and potassium 
halides are not very striking in view of the facts 
that the spectra were recorded in different labo- 
ratories and are based on measurements at two 
different wave-lengths in which the dispersions 
obtainable, in terms of energy units, differ 
greatly. With reservations, it can therefore be 
said that the pattern in Fig. 7 confirms the in- 
verse proportionality with the squares of the 
lattice distances, as required by Kronig’s 
hypothesis. 

On the other hand, Brewington’s data for the 
absorption structures accompanying the K-ab- 


sorption edge for Cl- in the various alkali 
chlorides show very significant variations from 
the relations predicted by Kronig’s hypothesis. It 
is conceivable that there may be some principal 
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difference between the secondary structure ac. 


companying absorption by anions as compared 
with cations, such as different selection rules for 
the electronic transitions.’ The fact, discussed jn 
the foregoing section, that the energy difference 
between presumably similar emission and ab. 
sorption transitions is not the same for Nat and 
Cl- transitions in NaCl, or for K+ and Cl- transj. 
tions in KCl, is in harmony with the supposition 
that electron transitions involving anions in ionic 
lattices obey different selection rules from those 
involving cations. Additional experimental data 
concerned with the secondary absorption struc. 
ture accompanying absorption edges of anions 
would be desirable. It is suggested that a study of 
the secondary absorption patterns accompanying 
the K-absorption edge for fluorine in the series of 
alkali fluorides, whose lattice constants extend 
over a wide range, would serve to make any 
differences very evident when compared with the 
present patterns for the sodium halides. 
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\ 


A radioactive gas possessing a half-life of 34.1+0.3 days has been produced by bombarding 
solid samples containing potassium, chlorine, calcium, or sulfur. The observed data require 
that this activity be assigned to A® in accord with the following reactions: 


Het 
S*+ He++A" +n 
Ca®+n—-A"+ Het. 


Alpha-bombardment of chlorine has failed to produce a measurable quantity of this activity by 
an (a, d) reaction. The detection of 4.72A Cl Ka x-rays by means of critical absorption measure- 
ments in gaseous compounds of chlorine and sulfur reveals that A® decays by K-electron 
capture. Cloud-chamber observations and aluminum absorption measurements indicate that 
K-electron capture accounts for more than 99.9 percent of the total number of disintegrations. 
The probability of the x-ray quantum being internally converted has been calculated to be 


0.96+0.03. No gamma-rays were observed. 


INTRODUCTION 


HE first radioactive argon to be reported! 
was produced by the bombardment of the 

gas itself with deuterons. This activity, which 
decayed with a half-life of 1.83 hours, was 
assigned to the isotope A‘'. The radiation emitted 
was found to consist of beta-rays? with a maxi- 
mum energy of 1.5 Mev and of a gamma-ray’ of 
1.37 Mev. More recently, a 1.9-sec. radioactive 
argon A** has resulted from the bombardment of 
chlorine with protons‘ and of sulfur with alpha- 


'A. H. Snell, Phys. Rev. 49, 555 (1936). 

*F. N. D. Kurie, J. R. Richardson, and H. C. Paxton, 
Phys. Rev. 49, 368 (1936). 

‘J. R. Richardson and F. N. D. Kurie, Phys. Rev. 50, 
999 (1936). 

‘*M. G. White, E. C. Creutz, L. A. Delsasso, and R. R. 
Wilson, Phys. Rev. 59, 63 (1941). 


particles.’ This isotope was found to emit posi- 
trons with a maximum energy of 4.38 Mev. 

The stable isotopes of argon have the mass 
numbers 36, 38, and 40. The existence of isotopes 
A*® and A** had never been proved although it is 
evident from the nuclear chart, shown in Fig. 1, 
that they might be produced by the bombard- 
ment of solid samples containing one of the 
neighboring elements. However, if an isotope of 
argon formed in this manner should decay en- 
tirely by K capture, then the resulting Ka x-rays 
of chlorine (4.72A) would be largely absorbed in 
the sample itself. It is well known that a sensitive 
means for detecting a soft radiation emitted by 
a radioactive gas is obtained by admitting the 


a 941) D. P. King and D. R. Elliott, Phys. Rev. 59, 108A 
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Fic. 1. Nuclear chart in the region of argon. 


gas directly into an ionization chamber. It 
seemed of interest, therefore, to perform experi- 
ments in which the radioactive argon is extracted 
from the bombarded sample and conducted into 
an ionization chamber. 

The 34.1-day radioactive gas which was found 
by this procedure has already been reported*® 
briefly and has been assigned to the isotope A*’. 
A thorough study of the radiation emitted by 
this isotope has now been made, and the results 
of this investigation are presented here, along 
with a more detailed description of the previous 
experiments leading to the isotopic assignment. 


TaBLE I. Relative yields from various reactions. 


Nuclear reaction 


K29(d, 
Cl*7(d, 2n)A” 

a)A® 
Cl87(d, 2m) 
Cl*7(p, 
S*(a, 2) A” 
Ca**(n, a) 

NaCl+Het 


*P. K. Weimer, J. D. Kurbatov, and M. L. Pool, Phys. 
Rev. 60, 469 (1941). 


ISOTOPIC ASSIGNMENT AND DECAY PERIOD 


Bombardments were made in the 42-inch 
cyclotron of The Ohio State University. A small 
quantity of the substance to be bombarded was 
placed in a copper target holder which was in- 
serted between the dees. A 2.40 mg/cm? alumi- 
num foil was placed in front of the sample to 
lessen contamination, but no attempt was made 
to seal the target to prevent the escape of 
the radioactive gas formed as a disintegration 
product. 

The gas was extracted from the bombarded 
sample by dissolving the sample, then bubbling 
air through the solution and allowing it to pass 
into an ionization chamber which had _ been 
previously evacuated. Suitable precautions were 
taken to prevent volatile radioactive substances 
other than rare gases from entering the chamber. 
The ionization chamber, which had a volume of 
about two liters, was connected to a Wulf bifilar 
electrometer. 

A summary of bombardments with the relative 
yield of each is given in Table I. Since the activity 
was obtained from a deuteron bombardment of 
potassium, the active isotope might be either A”, 
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A®*, or possibly A**. As the 34.1-day activity was 
obtained by bombardment of chlorine with 
tons, the activity must therefore be assigned 
to the isotope A*’. The remaining bombardments 
listed in Table I are all consistent with this 
assignment. Further evidence that the 34.1-day 
activity is not owing to A** was given bya neutron 
bombardment of potassium which produced only 
the 1.83-hour A“. 
The decay curves for each bombardment are 
shown in Fig. 2. The half-life was observed to be 


-34.14+0.3 days with no indication that any 


other period was produced by any of the bom- 
bardments. 


PROOF OF K-ELECTRON CAPTURE 
Critical Absorption Measurements 


If A®’ decays by K-electron capture, the char- 
acteristic 4.72A chlorine Ka x-rays would be 
emitted. This wave-length may be identified by 
the fact that it will be very strongly absorbed in 
sulfur but only slightly absorbed in chlorine. 
Because of the high absorption coefficients in- 
volved, it was considered most convenient to 
use the gases hydrogen sulfide and freon (CCI2F 2) 
as absorbers. 

For this purpose a special ionization chamber, 
shown in Fig. 3, was constructed. The radio- 
active argon was confined at atmospheric pres- 
sure in the small chamber marked A. A thin 
aluminum window of 2.4 mg/cm? thickness sepa- 
rated the active gas from the absorbing gas 
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Fic. 2. Decay curves for argon A®’. 
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Fic. 3. Gas absorption ionization chamber, showin 
movable screen S, inner cylinder C, and position of sma 
gas chamber A. 


which was placed at atmospheric pressure in the 
ionization chamber. The inner cylinder C and 
the movable copper screen S were kept at 135 
volts negative with respect to chamber A and 
the central collecting rod. Hence, only the ions 
formed below the screen had their charge col- 
lected by the central rod and measured by the 
electrometer. Thus by adjusting the distance of 
the screen from the foil, the thickness of the 
absorbing gas could be varied. 

The advantage of this method is that it is 
unnecessary to vary the pressure of the absorbing 
gas as is done in the usual methods of measuring 


gas absorption. Therefore, a much thinner foil . 


could be used for separating the active gas from 
the absorbing gas than could be used otherwise. 
This advantage is very great since an aluminum 
foil with a thickness of only 2.40 mg/cm? will 
reduce the intensity of a 4.72A x-ray to 8 percent 
of its initial intensity. 

The collecting rod of the gas absorption 
chamber was connected to a sensitive Wulf 
unifilar electrometer. A 12-hour deuteron bom- 
bardment of KCI was found necessary to produce 
sufficient argon to give a measurable activity 
through the aluminum window. Even then, the 
activity was weak and amounted to only 1.5 
times background when the screen was closest 
to the foil. 
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In Fig. 4, the ionization current obtained when 


chamber A of Fig. 3 was filled with argon A® is 
plotted against the corresponding screen posi- 
tion. The more rapid initial drop of the ionization 
current as the screen was lowered in hydrogen 
sulfide (curve C) than was observed in the case 
of freon (curve B) indicates that the absorption 
coefficient in hydrogen sulfide is the greater in 
spite of the fact that its density is less. However, 
when an additional 1.7 mg/cm? (0.25 mil) 
aluminum foil was placed over the aluminum 
window of chamber A, curves D and E show 
identical absorption in freon and in hydrogen 


sulfide. 


The decrease in ionization current as the 
screen is lowered is not entirely owing to absorp- 
tion in the gas since at the same time the effective 
part of the ionization chamber is being decreased 
in volume and is being moved a greater distance 
from the radioactive source. The effect of these 

. factors alone may be studied by filling chamber A 
with a radioactive gas emitting a much harder 
radiation. A radioactive gas which has been 
found to be convenient for this purpose is the 
1.44-day krypton Kr7** which emits beta-rays’ 
with maximum energy of 0.4 Mev. With chamber 


1.0 A. Nitrogen 
B. Freon (C Ci,F,) 
C. Hydrogen Sulfide 
09 D. Freon + .25 Mil Aluminum 
s 0.7 
§ 
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_ Fic. 4. Gas absorption curves for argon A*’, showing 
ionization current as a function of the position of the 
movable screen S. 


A containing krypton, the ionization currents 
were measured when the screen was lowered in 
freon, hydrogen sulfide, and in air, and curves 


7E. C. Creutz, L. A. Delsasso, R. B. Sutton, M. G. 
White, and W. H. Barkas, Phys. Rev. 58, 481 (1940). 


were plotted in a similar manner to those shown 
for argon in Fig. 4. In the case of krypton, the 
ionization current decreased less rapidly as the 
screen was lowered, and the curves for freon, 
hydrogen sulfide, and air were identical within 
the limits of experimental error. 

Since a 1.7 mg/em?* aluminum foil should re. 
duce an x-ray beam of 4.72A to 16 percent of its 
initial intensity, it is apparent from Fig. 4 that 
in addition to the x-rays, there is a harder com. 
ponent present in the radiation which passes 
through the aluminum window. Hence, by sub. 
tracting the ordinates of curves D from B, and 
E from C, the approximate intensity of the x-rays 
alone may be calculated at each screen setting. 
After making a correction for the dependence of 
the sensitivity of the ionization chamber upon 
the screen position, the experimental adsorption 
curves for the x-rays in freon and in hydrogen 
sulfide were plotted on semilog graph paper. 
These. curves were in satisfactory agreement 
with the calculated absorption curves for 4.724 
x-rays in these gases, with the absorption in 
hydrogen sulfide greater than that in freon by 
about the expected amount. Only x-rays lying 
between the critical absorption limits 4.38A and 
5.008A will possess absorption curves in these 
gases which compare even approximately with 
the observed absorption curves. If the soft 
component of the radiation passing through the 
foil were owing to very low energy electrons, then 
a larger absorption coefficient should have been 
observed in the more dense freon. 

It is therefore concluded from the critical 
absorption measurements that the characteristic 
Ka x-rays of chlorine have been identified and 
that K capture in radioactive argon A” is 
confirmed. 


Aluminum Absorption Measurements 


Aluminum absorption measurements on the 
radiation emitted by the argon were taken after 
the gas absorption experiments indicated the 
presence of an additional harder component 
passing through the foil. The absorption curve 
obtained from these measurements is shown in 
Fig. 5. The data for the points indicated by open 
circles were taken with the gas absorption ioniza- 
tion chamber by placing additional aluminum 
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foils below the window. The data for the solid 
circles were obtained by mounting chamber A 
above a more sensitive ionization chamber filled 
with freon at 20-lb. pressure. The resulting curve 
shows clearly that there are two components 


© Pressure tomzeten Chember 


Specie! 
Subtrecton of Curves & ong 
Curve C Cotewleted for X-Reys 4.724 


Mg Cm Aluminum 


Fic. 5. Aluminum absorption curve for argon A*’. 
Section of curve labelled A shows effect of x-ray component, 
while B is owing to the harder component of range 95 
mg/cm? in aluminum. Curve C represents the absorption 
in aluminum of 4.72A Cl Kae x-rays. 


present in the radiation passing through the 
window. The points obtained after subtracting 
the harder component give a curve consistent 
with the calculated absorption curve for 4.72A 
x-rays in aluminum. 

The range in aluminum of the harder com- 
ponent was observed to be 95 mg/cm’. Siace the 
intensity of the radiation was weak, it was not 
possible to determine whether this component 
was composed of positrons or electrons. How- 
ever, if we assume that the component does con- 
sist of one or the other, then the range in alumi- 
num is found to correspond® to an energy of 
0.33 Mev. Because of the high specific ionization 
of the beta-rays and the high internal conversion 
coefficient of the x-rays as well as their large 
absorption in the aluminum window, it follows 
that the harder component must represent a 
very small fraction of the total disintegrations 
occurring inside chamber A. The results of 
Table II, discussed in a later section, can be 
used to calculate this fraction, and it is found 
that regardless of whether the harder component 


*P. Curie, Radioactivité (Hermann and Company, Paris, 
1935), p. 294. 


is owing to positrons from A® or to a beta-emit- 
ting impurity, the component cannot account for 
more than 0.1 percent of the total disintegrations. 


Cloud-Chamber Observations 


A small quantity of radioactive argon mixed 
with air and hydrogen was admitted into a cloud 
chamber. A typical photograph, shown in Fig. 6, 
reveals many spots and clusters of drops but no 
electron tracks. Similar spots appearing in cloud- 
chamber photographs have been interpreted by 
others® as owing to short range photoelectrons 
ejected by soft x-rays. Such spots or clusters prob- 
ably represent single disintegrations. Clusters 
may be formed instead of spots when the time 
interval between the disintegration and the ex- 
pansion is sufficient for diffusion of the ions. 

When A*® decays by K-electron capture, an 
x-ray photon of 2600 electron volts is emitted 
by the atom provided the photon is not internally 
converted in an outer electron shell. In case con- 


Fic. 6. Cloud-chamber photograph of argon A*® ob- 
tained by admitting the gas into expansion chamber. 
Spots and clusters of drops are due to Auger electrons and 
short range photoelectrons. 


version does take place, an Auger electron will 
be ejected from the atom with an energy ranging 
from 2400 to 2600 volts. An electron of this 


*E. J. Williams and E. Pickup, Nature 141, 199 (1938). 
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energy has a range in air of about 0.04 cm. Thus 
a spot or cluster may be produced by an Auger 
electron emitted by the disintegrating atom 
‘itself as well as by a photoelectron ejected from 
a neighboring atom by the x-ray. 


INTERNAL CONVERSION COEFFICIENT OF 
X-RAYS FROM K CAPTURE 


When radioactive argon mixed with air is ad- 
mitted directly into an ionization chamber, each 
disintegration by K capture releases 2600 ev of 
energy regardless of whether the x-ray is in- 
ternally converted or not. However, when the 
gas is enclosed in chamber A and the ionization 
is owing to the radiation penetrating the foil, only 
those atoms in which the x-ray quantum is not 
internally converted are contributing to the 
ionization current. The Auger electrons, having 
an energy of less than 2600 ev, are entirely too 
weak to penetrate the 2.4 mg/cm? aluminum 
window. It is then apparent that a possibility 
exists for calculating the internal conversion 
coefficient for the x-rays by observing the ioniza- 
tion current produced through the foil by the 
argon, and then admitting the same gas directly 
into an ionization chamber. When the argon is 
in chamber A it is, of course, necessary to correct 
for the absorption of the x-rays in the aluminum 
window and for the less advantageous geometry. 
In order to make a better estimate of the geo- 
metrical factor, the measurements for argon have 
been compared with those for krypton Kr *. 

In Table II are recorded the activities of argon 
and krypton each measured first in chamber A 
and then directly inside an ionization chamber. 
The value of R represents the factor by which 
the ionization current is multiplied when the 
gas is placed inside the ionization chamber. In 
the case of krypton, where the radiation consists 
of 0.4-Mev positrons, the absorption in the 
aluminum window may be considered negligible, 
and the value of R is determined largely by 
geometrical conditions. However, for argon R is 
influenced also by the absorption of the x-rays in 
the aluminum window and by the internal con- 
version of the x-rays. The probability that the 


TABLE II. Relative activities of A and Kr. 


Argon Krypton 
T Half-life (days) 34.1 1.44 
t Bombardment time (hours) 12 2 

I, Activity of gas in 

chamber A 0.060 1.13 
Iz Activity of same gas in 

ionization chamber 415 38.6 
R Ratio of J, to 1; 6910 34.2 


x-ray quantum will not be converted equals 
approximately (34.2) (0.60) /(6910) (0.08) = 0.037. 
The quantity (0.60) appears because only 60 
percent of the radiation observed through the 
window of chamber A is composed of x-rays while 
the factor (0.08) corrects for the x-ray absorption 
in the aluminum window. Therefore, the value 
of the internal conversion coefficient which repre. 
sents the probability of conversion is found to be 
0.96+0.03. 

The internal conversion coefficient of the Ke 
x-rays of chlorine has been previously deter- 
mined by means of fluorescent yield experiments 
in which the K electron is removed by irradiating 
the atom with x-rays. Although it may be 
questionable whether the internal conversion 
coefficient is the same when the K electron is 
captured by the nucleus as when the electron is 
ejected photoelectrically from the atom, a com- 
parison of the values calculated by the two 
methods is of interest. An average value of the 
coefficients obtained by several observers!® using 
the fluorescent yield method is 0.92, which is 
somewhat lower than the coefficient 0.96+0.03 
arrived at in this paper. A value also greater than 
0.92 was calculated from theoretical considera- 
tions by Burhop." 

It is a pleasure to acknowledge the support 
received from the University Development Fund 
and from Mr. Julius F. Stdne. Valuable assistance 
received through the graduate school is also 


acknowledged. 


10 A, H. Compton and S. K. Allison, X. in Theory and 
Van Nostrand Company, Inc., New York, 
1935), p. A 
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Probability equations are developed for the purpose of analyzing recrystallization data to 
determine what mechanisms may operate during recrystallization processes. The equations 
allow for inexact variations in orientation changes provided the variations are small. 
P(E) = 1.43e' gives the probability that a nucleation will result in an orientation falling within 
a radians of a specified one A on a chance basis. Other equations are given for the probability 
of occurrence of nearly common orientations for various types of groupings. Values of prob- 


abilities are given in tables. 


INTRODUCTION 


EW crystals are produced whenever a single 
crystal is deformed a small amount and 
heat-treated above its recrystallization tempera- 
ture. In general, the changes in orientation that 
occur cannot be described in terms of simple 
mechanisms such as twinning or as the rotation 
of the lattice about a single axis. Furthermore, 
if a mechanism is postulated that involves a 
large number of transformations, it may be of 
little use in predicting changes in the crystal 
lattice unless appropriate selection rules are also 
introduced to limit the number of final orienta- 
tions. Such rules would depend upon experi- 
mental conditions. 

In an attempt to overcome some of these dif- 
ficulties, the writer has found it useful to analyze 
recrystallization data with the aid of probability 
equations. These equations are developed in the 
present paper. The first relationship, Eq. (1), 
gives the probability that a single nucleation will 
produce an orientation falling within a radians 
of a specified one A assuming all possible orien- 
tations as equally likely. Any nucleation that 
satisfies this condition is called event E and its 
probability, P(Z). The reciprocal of P(£) will 
be called NV. It is the maximum number of events 
similar to E. We might say that there are N ways 
for an orientation to occur and if each way is 
equally likely, then the probability that it is E 
is 1/N. Call each of the N configurations a state 
(a classified orientation). These states, instead of 
orientations, may be used in considering’ mech- 
anisms of recrystallization. The number of states 
allowed for a given mechanism, however, must 
be considerably less than N if the mechanism is 


to have much meaning. In this connection, an 
illustration is given for a mechanism that would 
produce high order twins. Finally, an application 
is made of P(E) to obtain Eqs. (2a) and (2b) 
for the probability of occurrence of new crystals 
having nearly common orientations. 


PROBABILITY P(E) OF OCCURRENCE OF AN 
ORIENTATION IN A GIVEN STATE 

Consider all orientations that agree within a 
radians of A, defined by three mutually per- 
pendicular poles A;, A2, As; of {100} planes, as 
belonging to the set £, i.e., in a state E. An 
orientation B belongs to E if no {100} pole of 
B departs from A by more than a radians. We 
ask, “What is the probability of occurrence of 


-event E if an orientation is selected at random 


from a group of random orientations?” If E is to 
occur, B must be taken from a subset of the 
group of random orientations having one pole 
within @ radians of one pole of A. Consider the 
subset for pole A; and call a selection from it 
event F. Complete a set F by taking all such 
orientations. From this set make further selec- 
tions, each having a pole within a radians of A, 
and thus obtain a set G. Finally, all orientations 
of this group that have a pole within a radians of 
A; may be taken to form a set H. 

The probability of E (event E) will be 
P(E) =P(F)P(G)P(H) where P(F) is the prob- 
ability of F, P(G) is the probability of G after F 
has occurred, and P(#) is the probability of H 
after F and G have occurred. »** 

‘Arne Fisher, Mathematical Theory of Probabilities, 
(The MacMillan York, 1922), p. 29. 

C.F ity and Its Engineering Uses 


Van Nostrand Company, Inc., New York, 1928), 
(a) p. 116; (b) pp. 19, 20. 
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/ In order to obtain P(E), the original group of arc equal to 2[(Ra)?—(R8)*}. Therefore, 4 
random orientations will be represented by a i! 
H plot of {100} poles on the surface of a hemisphere 2| 2k(a—8)2[(Ra)?—(RB)*}'Rdg 
i of radius R. Let k be the number of poles per P(G)= 0 
unit area, thus giving a total of 27R*k poles and 
2xR*k/3 orientations. Computing the area near 2 2k(a—f)x/2 R Rdg 
A,, we find that there are 2x(1—cos a)R*k poles ( 
within @ radians of A;. Each pole belongs to a 2a(1—4/3x) =1.15a. le 
! unique orientation. Consequently, there are + P a 
2x(1—cos a)Rtk orientations that belong to the tribured within ating 
set F. Expressing probability in the usual way, ions 2eR in a rather complicated way. Con. 
\) we get at once for low values of a: sequently, the integration required to obtain 
2x(1—cos a) would be difficult to carry out. However, 
i P(F)= = 3(1—cos a) = ja’, 
2eR*k/3 
yt the error being only 74 of one percent for a 
| equal to 10°. 
| Figure 1 is a stereographic projection of the 
y plot of set F. The poles have a density & in the or 
large circle about A. The distribution of poles in he 
the band of width 2a@R containing Az and A; is th 
I assumed to be uniform along the band, but the a 
i distribution across the band must be found. an 
I Divide the large circle about A; into rings of 
i width Rdg. Consider one of radius Rg. It contains 
i approximately 2rR8RdBk poles, the error being TI 
i 3 of one percent for a equal to 10°. Twice this als 
number of associated poles are spread uniformly eit 
i along 180° arcs that fall within a band of width 
: 28R (see Fig. 1). These arcs produce a uniform 
i distribution of the poles. Considering every Fic, 1. Stereographic projection of the plot of set F together | 
(t value of 8, it may be seen that the distribution with diagrams used in calculating P(G). al 
i across the band, to a good approximation, is ; | po 
: linear, varying from zero for 8 equal to a toa WS ©" Say definitely that in 
f maximum at the center of the band. It may be ro? R? of. 
represented by a(a—) and must satisfy the fol- (2a (H) <1, we 
sta 


! lowing equations, the integration being carried 


out over } the band along a 90° arc: or P(H) is between +/4 and 1. Substituting the 


values for P(F), P(G), and the limits of P(H)in | eit 


} 

rot R= 2 Rag. the equation for P(Z), we find that val 
1.35a*< P(E) <1.72a°. 
Rr/2 2aR gives the approximate area wherein Although this range is sufficiently accurate for ing 
f one-half the associated poles fall. Consequently, _ most applications, a more definite value of P(E) wa 
' a=2k. may be obtained as follows: of 
| We must now obtain all the poles in the circle write P(E) =P(F)P(0). ~ 


about A: in order to calculate P(G). Consider an 
arc where the density is 2k(a—8). To a good The set J is to be obtained from F in sucha _ for 


approximation the circle intercepts a length of way that all poles fall within the circles about to 


i 
it ° 
| 


| 


ire dis. 
dimen- 
Con- 
obtain 
wever, 


A; and As3. Since a lengthy derivation is needed 
in this case, only the result will-be given, namely : 


32 
P(I) —V})a=0.95a 


(as before, approximations have been made for 
length of arc). Substituting the values for P(F) 
and we get: 


P(E)= 1 4303, (1) 


As an illustration, Eq. (1) gives P(E) = 1/1140 
for a equal to 5°. In this case, the maximum 
number of states (JV) is 1140. 


PROBABILITY OF OCCURRENCE OF AN 
ORIENTATION IN ONE OF SEVERAL 
GIVEN STATES 


Case I: Two states, K and L, in which all 
orientations are uniquely classified. (We may 
have overlapping circles in this case, but one of 
the three circles for one state must not overlap 
a circle of the second state). The probability that 
an orientation will occur in state K by chance is: 


P(K) =1.43a'. 


The probability that it will occur in state L is 
also 1.43a*. The probability that it will occur 
either in state K or L is: 


P(K or L) = P(K)+P(L) =2X1.43a'. 


Case II: Two states, S and T, overlap so that 
a unique classification of all orientations is im- 
possible and some have to be described as being 
in states S and T. In this case, we should think 
of each as a state only when it is alone; together, 
we should think of them as being less than two 
states. 

The probability that an orientation will occur 
either in state S or T by chance is less than the 
value obtained for Case I. We note that if S 
and T are allowed to merge, the probability 
drops from 2 X 1.43a* toward 1.43a', correspond- 
ing to a transition from two states (Case I) to- 
ward one state. A correction owing to the amount 
of overlapping area might be made for the circle 
having the smallest amount of overlap, but the 
non-uniform distribution of poles within a circle 


_ for a state would make this correction difficult 


to determine. Nevertheless, there would be a 
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number f such that P(S or T)=f1.43a*, and 
1<f<2. 

Case II]: When data for standard twin groups 
through three generations are plotted on a 
stereographic net,* it is found that no over- 
lapping occurs among the 52 orientations until 
a becomes equal to about 23°. The probability 
that a new crystal will occur in one of the 52 
states, when all N states are equally likely, is: 


2.59r\3 
P=52/N=52P(E) =52X1.43 (=) 


= 0.00619 = 1/161. 


When 5° is taken for a, the number of states 
N is 1140, but the maximum number of states 
allowed on the basis of a mechanism that pro- 
duces only three orders of twins is 52. Because 
of overlap, there is some degeneracy among the 


states and their effective number is reduced as. 


pointed out in Case II. Consequently, 


This result means that recrystallization with- 
out any mechanism or preferment should not 
result in a 5° type fit more than once in every 22 
times if we restrict a fit to one of the 52 first-, 
second-, and third-order twins. A mechanism of 
recrystallization of this degree of complexity 
should not be hard to test experimentally.‘ 


PROBABILITY OF OCCURRENCE OF ORIENTA- 
TIONS FOR VARIOUS. TYPES OF GROUPINGS 


Consider a group of random orientations exist- 
ing in N states, Hi, He, H3---Hw (N=1/P(E)). 
Select n orientations (m< N) in an independent 
manner and suppose that these fall into / states 
in such a way that m, belong in one state, mz in 
another state, etc., up to m, in a kth state (m2 2), 
and one each belongs in the remaining /— states. 
Suppose further that some of the & states con- 
tain an identical number of orientations. Desig- 
nate the number of equal size groups by 11, 72, 
etc. We ask, ‘What is the probability of occur- 


* Data on four generations of twins have been obtained 
at the Pittsfield Works’ Laboratory and are now being 
for publication. 
‘In this connection, considerable data have already 
been accumulated at the Pittsfield Works’ Laboratory of 
the General Electric Company. 


52 1 
P<— or —. 
| 1140 22 
ng the 
(H) in | 
ate for 
P(E) | 
uhe | 
about 


218 CECIL G. DUNN 


rence of such a group?” The group in this case The number of combinations, however, is jp. 


may be represented as follows: creased if all the m’s are not the same. Suppose = 
f _ one of the m’s occurs 7; times, another r2 tim 
i m(1) (k) 1 1 €s, 


ia tations in a particular state. We do not, however, The number of combinations fending, they 
| 


specify the states in advance. As an example, fore, to a 

consider the case of throwing three dice to obtain 

if a pair. In one throw we use three orientations ¢ 7" 

(n) out of six (VV), but to get a pair one orienta- ype 


tion must occur twice. Therefore, / equals 2, _— 
k equals 1, and m, equals 2. Figure 2 gives an (N—Rk)!k! ry !re!--- 
illustration of a more general orientation group 


or 


We may also write combinations for the /—, 
states that are occupied by only one orientation 
each. Since k states of the N have been used in 
forming the group m,m2---m,, there are N—} 
remaining to select from to get the /— singly- 
occupied states. The number of combinations is: 


(N—k)! (N—k)! 
(N—J)!(—k)! 
Any one of the combinations containing n orien- 


tations could have been selected in a number of 
unique ways, namely, 


Therefore, the total number of ways to form a 
specified type of group is: 


| : Fic. 2. Stereographic diagram showing states 1, 2, 3, 4, N! (N—k)! - mt 
{ and 5 containing 3, 2, 2, 1, and 1 orientations, respectively. : = 
a is 5°. N21000. 

k 
; existing in five states with m, 1, k; mi, m2, ms; —_ Since the total number of ways to make n ~ 
" r, and re equal to 9, 5, 3; 3, 2, 2; 2 and 0, re-_ selections is N*, the probability of getting the : 
i spectively. specified type of group is: 2 
i To obtain the probability of occurrence for the ; 
: general case, consider first any group of k states i 2 
| of the following type: Nin! ‘ 
m m m (2a) 
The number of combinations of N states taken | When /=n, this equation reduces to: “ 
: k at a time for a group where NI | 
mM, =m (N—n)!N* P, 


(N—k)!k! 


*In superscripts and subscripts m; is written as m/(1); 
mz as m(2); mz as m(3); etc. 


an expression giving the probability that all m 
orientations occur in different states. 


i 
| 
@' QO} 
| O)s Qs Os 
C.)3 
n! 
— 
| 


(2a) 


all 
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TaBLeE I. Probability of occurrence of orientations for various types of grouping. 


of 
b a = a ms Probability from Eqs. (2a) and (2b) Probability 
N! 
0 Pe 
(N—n)!N* 
Pe (N—n)!n! P® n(n—1) 
1 2 
2! (N—n+1)'(n—2)! 2 N—n+1 
P®  (N=n)!n! Ps (n—2)(n—3) 
2 2 2 
212! (N—n+2)!(n—4)!2! N—n+2 
1 P® (N=—n)!n! 1 Ps (n—2) 
1 3 P3=— 
3! (N—n+2)!(n—3)! 3 N—n+2 
3 2 2 2 22> 2=— 
21212! (N—n+3)!(n—6)!3! 8 N—n+3 
P® (N=—n)!n! Ps (n—3)(n—4) 
2 2 3 
312! (N—n+3)\(n—5)! 2 N—n+3 
P® (N=n) In! Pas (n—5)(n—6)(n—7) 
(2!)4 (N—n+4)\(n—8) 32 N—n+4 
1 Pp (N—n)!n! 1 P; (n—3) 
1 4 
4! (N—n+3)\(n—4)! 
P® (N=—n)!n! 1 Py (n—4) 
1 
5! (N—n+4)!(n—5)! 5 N—n+4 
At least two in one state 1-—P° 1—P°® 
k 
TABLE II. Values of Pa(i)m(2)-+ for N=1000, 
a of n=20 n=10 n=5 
0 0.8259284  0.8259284 0.9558606 0.9558606 0.9900350  0.9900350 
0.1599657  0.9858941 0.0434044  0.9992650 0.0099401 0.9999751 
0.0124616  0.9983557 0.0006125  0.9998775 0.0000149 —0,9999900 
0.0009773 0.9993330 0.0001166 0.9999941 0.0000099 0.9999999 
0.0005070  0.9998400 0.0000030  0.9999971 — 
0.0001352 0.9999752 0.0000024 0.9999995 9.9x107 
na 2 0.0000117 0.9999869 4.5X 10 
0.0000042  0.9999911 0.0000003  0.9999998 5.0107 
1.4 10-* 3.21071 9.9X 10-8 
At least two in one state 0.1740716 0.0441394 0.0099650 


If we substitute P® in Eq. (2a), we get: 


m(2)+++m(k) 


P°(N—n)!n! 


Successive applications of Eq. (2b) give the 


probabilities listed in Table I. 


If a small value of n is taken, it can be shown 


readily for any N whatsoever that the sum of 
the probabilities for all possible groups is one. 
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For instance, 


when is 5 and NS5. 

The values of N!, etc., in the expression for P® 
may be obtained by using Stirling’s formula. If 
n is small, however, it is easy to compute P* 
directly with the aid of either logarithms or an 
electrical computing machine. 


N. FERGUSON, JR. 


Table II gives some numerical results op 
probabilities for one value of N and three 
values of n. 
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The Photoconductivity of Sodium Chloride in the Far Ultraviolet 


J. N. Fercuson, Jr.* 
Cornell University, Ithaca, New York 


(Received August 7, 1944) 


Measurements have been made on the wave-length dependence of the photoconductivity of 
an uncolored sodium chloride crystal in the far ultraviolet region of 900 to 1350 angstroms in 
order to test the theories of Peierls, Frenkel, and Mott regarding the existence of discrete 
energy levels, called exciton levels, lying just below the conduction bands. A region was found 
between 1300 and 1350 angstroms in which there was large absorption but very small if any 
photoconductivity or external photoelectric effect indicating the existence of exciton states. 
At wave-lengths shorter than 1300 angstroms both the photoconductivity and the photo- 
electric effect began to set in and these showed a doublet structure. The photoconductivity 

_ rises abruptly at the second absorption peak in the absorption spectrum and shows structure 
which locates energy levels above the ground state. The experiments of Podubnij on the 
iodides are a contradiction to this experiment inasmuch as the alkali halides should reasonably 


be expected to behave similarly. 


INTRODUCTION 


EIERLS,' Frenkel,? and Mott* have intro- 
duced the concept of non-photoconductive 
absorption in crystalline insulators because of the 
possible existence of discrete energy levels lying 
just below the conduction bands. The withdrawal 
of an electron by light of quantum energy corre- 
sponding to the distance AB, Fig. 1, leaves a hole 
in the filled band below A. This electron, termed 
an “‘exciton’’ by Frenkel, is then thought to be in 
the electric field of its positive hole, giving rise to 
discrete energy levels. Absorption of light of 
quantum energy corresponding to the distance 
AE would be a photoconductive absorption. A 
test for the existence of these exciton states would 
be to compare the dependence of photoconduc- 
* Now at Corning Glass Works, Corning, New York. 
1R. Peierls, Ann. d. Physik 13, 905 (19. 2). 


4 Frenkel, Physik. Zeits. Sowjetunion 9, 158 (1936). 
. F. Mott, Trans. Faraday 34, 500 (1938). 


tivity and absorption on wave-length for an 
alkali halide crystal. Experimental data on the 
photoconductivity of pure alkali halides are very 
meager. To the author’s knowledge the only 
previous work is by Podubnij‘ in which potassium 
and sodium iodide were illuminated over a wave- 
length range of 2500 to 1900 angstroms. The data 
showed photoconductive currents which depend 
in the same way on wave-length as does the 
absorption, and Podubnij concluded that all the 
absorbed energy is used to produce photocon- 
ductivity. There is no evidence that any care was 
taken to separate the primary photoconductive 
currents from secondary currents which Gudden 
and Pohl® have shown must be done in order to 
obtain interpretable results. 
* Podubnij, J. Exper. and Theo: Phys. 8, 410 (1939). 


5 Hughes and DuBridge, Photoelectric Phenomena (Mc- 
Graw-Hill Book Company, Inc., New York, 1931), p. 291. 


é 

| 

| 

| 

| 
ph 
tio 

dus 
sid 
i} the 
eas 
i AI 

] 
I by 

i 

| 


Its on 


ciation 
reham 
lebted 
Works 


1944 


PHOTOCONDUCTIVITY OF SODIUM CHLORIDE 221 


VACUUM 

CONOUCTION 
E Levers 

ExcitTON 

——p 

> 

O 

Zz 

uw 

10 

A 
Grouno sTATe(s) 


Fic. 1. Energy level scheme. 


In view of the meager information about the 
photoconductivity of pure alkali halides, in addi- 
tion to the desirability of obtaining evidence to 
confirm the existence of the exciton states intro- 
duced by Peierls, Frenkel, and Mott, it was con- 
sidered worth while to attempt measurements of 
the photocondictivity of an easily produced and 
easily handled salt such as sodium chloride. 


APPARATUS AND EXPERIMENTAL PROCEDURE 


Illumination of the crystal was accomplished 
by mounting the crystal, as shown in Fig. 2, ina 


W.E 96475 


Leavs EvecTROmETER 


To Barreries 


Switch Box To CHance 
a GrRio CONNECTION From 
TO Caysrac. 


vacuum ultraviolet normal incidence spectro- 
graph. The grating could be rotated to give 
various wave-lengths on the crystal. The source 
of illumination was a capillary air discharge tube. 
An electrically controlled shutter was placed at F 
between the exit slit and the grating so that the 
crystal could be illuminated for short fixed 
periods of time. 

In order to measure relative intensities of the 
spectral lines used, a windowless silver cathode 
photoelectric cell was placed at E on the side of 
the crystal holder away from the grating. The 
crystal could be moved into or out of the beam 
from outside the vacuum. 

The electric circuit for measuring the quantity 
of charge released by light flashes and for 
measurements of currents in the photoelectric 
cell was the Barth type electrometer tube circuit 


whose sensitivity was 310,000 millimeters per 


volt drop across the input resistor. 

Thin single crystals of sodium chloride were 
produced by cleaving artificial crystals of high 
purity.* In order to decrease the possibility of 
surface conduction being measured as volume 
conductivity, the electrodes were behind the 
crystal and the crystal was mounted between 
them and the exit slit of the spectrograph. The 
technique developed by Gudden and Pohl’ for 
elimination of secondary effects of polarization 
and decrease of resistance owing to release of elec- 
trons in the crystal was followed. This consisted 


* Harshaw Chemical Company, Cleveland, Ohio. 
7 Hughes and DuBridge, reference 5, p. 289. 
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Fic, 2. Schematic diagram of apparatus. 
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in measuring all crystal photo-effects ballistically 
with low light intensity flashes and with a fairly 
low applied electrical field intensity in the crystal. 
The procedure was to set a spectral line on the 
slits by observing the photo-cell currents with the 
crystal out of the beam. After the intensity was 
observed to be steady within ten or fifteen 
percent, and it was usually much steadier than 
this, the shutter was closed; the crystal slid into 
position; and the electrometer tube grid con- 
nected to the crystal. The shutter was then 
opened for a fixed period of about one second by a 
timing device. The ballistic reading thus obtained 
was found to be made up of both external photo- 
electric currents as well as photoconductivity 
currents. These were separated in the following 
manner: The reading just described was made 
with a retarding potential, small compared to the 
potential drop across the crystal, between the 


E 


| 20 wioa | 
126 134 142 


10.2 
ENERGY [ELECTRON VOLTS) OF THE INCIDENT UGHT QUANTA 


Fic. 3. Plot of photoconductivity current and photo- 
electric effect current in arbitrary units per unit absorbed 
energy as a function of the energy in electron volts of the 
incident light quanta. 


brass shield, Fig. 2, and the crystal electrodes. 
The value of this potential was adjusted so as to 
be sure that all photoelectrons ejected from the 
crystal were collected by the positive crystal 
electrode. The reading i, is the photoelectric 


effect of the crystal 7,, plus the photoconductivity 


i-. To obtain measurement of the photoelectric 
effect alone, a small potential was applied across 
the crystal in a direction opposite to the large 


potential in the above measurement. This smal] 


potential was to insure that photoconductivity 
would not contribute appreciably to the current. 
A large accelerating potential was next applied, 
and a ballistic reading was made as before. This 
is the current 7,. Saturation of the current thus 
measured with accelerating voltage indicated 
that most of the photoelectrons were being meas. 
ured. No measurable deflections were obtained at 
any wave-length if the crystal were removed 
from the crystal holder and normal operating 
conditions established. 

Sources of error in this experiment are the 
following: (1) Unsteadiness of the light source. 
The light source gave two sources of unsteadiness;: 
one owing to the nature of the method of excita- 
tion of the discharge ; the other tonon-uniform rate 
of leak of air into the tube. The discharge was re- 
markably steady, however, and the photoelectric 
cell showed that over a period of half an hour the 
average intensity taken over the period of the 
galvanometer system never varied by more than 
25 percent and was usually 10 percent or less. 
(2) Non-monochromatism of the light flashed on 
the crystal. This was owing to three causes: first 
was the necessity of a broad slit width of ap- 
proximately 0.5 mm in order to furnish enough 
light to get substantial galvanometer deflections. 
The second was the scattering of light into the 
slit. This was of course minimized in the custom- 
ary ways by baffles and black walls. The third 
cause was the presence of overlapping orders of 
the spectrum. It is thought that the second-order 
intensity is small compared to the first-order 


intensity. Some evidence of this is that the photo- ° 


cell showed less than one-third as much current 
for the line 1085A in the second order as it did in 
the first order. (3) Unsteadiness of the electrome- 
ter tube. The proximity of a discharge tube to the 
electrometer tube brought about magnetic and 
electrical pick-up. These were so minimized that 
during the final runs reported here the tube was 
steady enough that readings could be made to 
plus or minus half-centimeter deflection. This 
was in general never greater than ten percent of 
the deflection. 
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PHOTOCONDUCTIVITY 


RESULTS AND DISCUSSION 


Overlapping data taken on separate days are 
plotted in Figs. 3 and 4. The photoconductivity 
current 7. Was obtained from the data by the 
relation 

Assuming that the yield of the photoelectric cell, 
je., the number of electrons released per unit 
incident energy, is a slowly varying function over 
the spectral region studied, one can take the 
photoelectric cell reading 7, to be proportional to 
the intensity of the light and compare intensities 
at different wave-lengths. This assumption was 
made and the yields in arbitrary electrical units 
per unit absorbed light energy were calculated, 
the crystal photoelectric effect yield being 7,/i,, 
and the photoconductivity yield being 7,/7,. It is 
noticed in Figs. 3 and 4 that there is a region of 
large absorption in which there is no measurable 
photoconductivity or photoelectric effect (even 
though the photo-cell showed considerable light 
intensity in this region), but that these both set 
in sharply and almost at the same wave-length. 
These facts are in agreement with the exciton 
theory; the absorbed energy in the region in 
which we have neither photoconductivity nor 
external photoelectric effect is used in the forma- 
tion of excitons. Curves for both the photocon- 
ductivity and photoelectric effect show doublet 
structure. Since the ground state of chlorine is a 
doublet of separation 0.11 electron volt, the first 
four peaks in the photoconductivity curve can be 
accounted for by assuming two narrow unfilled 
energy bands at 9.53+0.02 and 10.06+0.02 
electron volts above the higher ef the doublet 
ground states. The other peaks and region of 
continuum are explained by locating a band be- 
ginning at about 10.5 electron volts and another 
at about 11.6 electron volts above the higher of 
the ground state levels. Since the data_were taken 
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at larger energy intervals in the shorter wave- 
length regions, location of the latter two levels is 
approximate and there may be other levels as 
well. In comparing the two curves of photoelectric 
effect versus energy of the incident light quanta 
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PHOTOEFFECT (ARB UNITS) /ABS ENERGY 


Fic. 4. The total photoelectric effect in arbitrary units 
per unit absorbed light energy as a function of the energy 
of the incident light quantum. 


and photoconductivity versus the energy of the 
incident light quanta, it is noticed that peaks in 
the former curve occur where troughs of the 
latter curve occur, and vice versa. This would be 
expected if one realizes that any electron, given 
energy equal to or more than a certain minimum 
amount, can be ejected from the crystal as a 
photoelectron and is so ejected unless the value of 
the energy given is the correct amount to place 
the electron in a photoconductivity band. In that 
case the electron is a photoconductivity electron 
which moves through the crystal in the general 
direction of the field and the photoelectric current 
is small for the corresponding wave-length. 
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Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


New Researches in Gyromagnetism! 


S. J. BARNETT 


U ene, of California, Los Angeles, California, and California 
Institute of Technology, Pasadena, California 


September 20, 1944 


T the Strasbourg Réunion sur le Magnétisme,* where it 

devolved upon me to make the report on gyromag- 
netic phenomena, the President, Professor P. Weiss, said 
to me: “‘Your work has now become classical, but it must 
be extended to more substances and to the same substances 
in different states; and for this purpose it is necessary to 
increase the precision. And you are the man that has to do 
the work.” This extensive program was, in fact, already 
underway, but it is urgent that many others should par- 
ticipate in it. It is difficult because the theoretically im- 
portant differences between the gyromagnetic ratios are 
relatively small and the errors which have to be eliminated 
are seriously large. This investigation constitutes the initial 
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part of the program, and consists of new and com 
independent determinations, by means of the effect known 
as magnetization by rotation, or the Barnett effect, of the 
gyromagnetic ratios of a number of ferromagnetic syb. 
stances. Although it is now certain that the work can be 
greatly improved after the end of the war, important con. 
firmations of earlier work have already been achieved, and 
discrepancies have been reduced. 

The general method of electromagnetic induction by 
which the effect involved was discovered in 1914 was used, 
but with numerous and important modifications. The 
sensitivity was greatly increased by using two galya. 
nometers, the principal galvanometer acting as a fluxmeter 
and the other being actuated by it with the help of two 
photoelectric cells or, in about one-quarter of the work, a 
Moll relay. In most of the work, including all of that with 
the photoelectric cells, deflections were closely proportional] 
to flux changes; and only this work is reported in the tables 
below. The magnetic rotors were identical with, or similar 
in construction to, those used by S. J. and L. J. H. Barnett 
in Washington, in magnetometer measurements on the 
same effect. The sources of error were reinvestigated and 
were largely eliminated, most of them by processes used 
in the earlier work. In all the new work the rotors were 
driven at exactly 60 r.p.sec. by small synchronous motors, 
In the greater part of the work the rotor axis was hori- 
zontal, and the rotor was driven from either end at will, 
It appears that only by driving from both ends can all of 
the errors be completely eliminated. In the latest work the 
axis was vertical, and the rotor was driven from above 
only. An important improvement consisted in making 
the sensitivity and gyromagnetic measurements simul- 
taneously. 


TABLE I. Gyromagnetic ratios (p) from new determinations (¢/m =1.759 e.m.u.). 


Rotor No. of sets pXe/m Rotor No. of sets pe Xe/m 
Electrolytic iron III 15 1.025 +0.012 Cobalt I 16 1.072 40,006 
Electrolytic iron IV 12 1.031 +0.011 
Norway iron 4 1.028 +0.011 Hopkinson's iron-nickel alloy 19 1.019 +0.016 
Mean for soft iron 31 1.028 +0.002 Permalloy 6 1.053 40.006 
Steel III 10 1.039 +0.008 Cobalt-iron 22 1.029 40.009 

Cobalt-nickel 13 1.080 +0.020 
Nickel I 17 1.054 +0.028 Heusler alloy 21 0.989 +0.006 


TABLE II. Gyromagnetic ratios (p) of fe 


‘erromagnetic substances. Values of pe/m from the three extensive investigations 
in the author's laboratories (¢/m =1.759 e.m.u.). 


Mean of II and 1.035 1.040 1.046 
mean I and III Both +0.003 +0.002 +0.004 


Material Soft Hopk. Perm- Cobalt- Cobalt- Heusler 
Investigation iron Steel Nickel k alloy alloy Cobalt iron nickel alloy 
I Barnett 1.049 1.047 1.031 1.016 1.054 1.070 1.067 1.068 1.011 
II Einstein- 1.032 1.038 1.051 1.051 1.023 1.046 1.085 1.025 1.076 - 
de Haas 
Ill Barnett 1.028 1.039 1.054 1.019 1.053 1.072 1.029 1.080 0.989 
Mean I and III Barnett 1.038 1.043 1.042 1.018 1.054 1.071 1.048 1.074 1.000 


+0.002 +0.000 +0.001 +0.019 +0.006 


1.017 1.050 1.078 1.036 1.075 =e 
+0.001 0.004 +0.007 +0.012 +0.001 
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The means of the gyromagnetic ratios obtained as 
indicated above are given in Table I. A set consisted of 48 
(in some cases 64) observations made in a symmetrical 
manner and on a strict time schedule. 

In Table II these values are compared with those ob- 
tained in Washington from an investigation of the same 
effect but by the magnetometer method, and with those 
(II) obtained in Pasadena from an elaborate investigation 
of the converse effect (Einstein-de Haas effect). For the 

of this table all values for soft iron are combined 
into one mean, likewise all values for steel, etc. 

For I the authors claimed only a mean error for all the 
rolors not greater than 2 percent. It has now long been 
clear, and the present investigation makes it still clearer, 
that the mean error for the individual rotors in the Washing- 
ton work is less than this. The values in IT are considered 
by the author the most precise values of the gyromagnetic 
ratios hitherto obtained, those for Permalloy and soft iron 
being estimated to be correct to one-half percent or less. 
In view of the difficulties involved the three investigations 
agree remarkably well. The large value for soft iron in I is 
owing to the inclusion of an excessive result for a certain 
rotor, electrolytic iron II, and is not confirmed by II, as 
electrolytic iron IV was cut from the same rod, and elec- 
trolytic iron III was cut from similar material. The large 
value for cobalt-iron in I is not confirmed by either II 
or III. 

For important papers on the theory, which, so far as it 
has yet been developed, is in agreement with such values 
as are given above. See Gorter and. Kahn? and C. J. 
Gorter ;* also reports by Gorter and Kronig and R. Forrer* 
and the discussions by Kramers and others.* See also 
R. Forrer.” 

at Mts 1942 Los Angeles mecting. See Proc. Am. Acad. Arts 
Sci. 75, No. 5 

*The Rapports of thi were printed in Ry a, 
International Institute of Intellectual © Cooperation 
Central Service for Scientific ee. but their xX distribution 
has hitherto been prevented by, the 

4C. J. Gorter, Phys. Rev. 60, 836 (1941). 


‘Gorter and Keents and R. Forrer, Réunion (1939); 


Forrer, Comptes 207, 1390 (1938). 


A Strong Infra-Red Radiation from Molecular 
Nitrogen in the Night Sky 


Jor: Stessins, A. E. WuitrorD, AND P. Swincs 


Carnegie Institution Mount Wilson Observatory, 
ifornia 


September 20, 1944 


OME four years ago, in the course of photoelectric 

measures of the colors of stars and nebulae, Stebbins 
and Whitford found an intense infra-red radiation from the 
night sky which is scores of, if not a hundred, times as 
strong as the ordinary persistent auroral line at 5577A. 
That the new radiation must come from the atmosphere 
was shown by its variation with zenith distance, by its 
usual progressive decrease through the night, and by its 
variability from night to night and from season to season. 
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Renewed interest was awakened on July 20, 1944 when 
the radiation was doubly strong and was varying rapidly, 
as much as 10 or 15 percent in five minutes. No marked 
magnetic disturbance occurred on that night. It was then 


suggested by Swings that the radiation was undoubtedly — 


molecular in origin and very likely could be traced to the 
(0,0) band of nitrogen at about 10,450A. 

With this suggestion as a guide the wave-length was 
determined with the photo-cell and various thicknesses of 
a glass filter with a steep gradient in the transmission 
curve across 10,000A. The filters were readily calibrated 
in the first-order solar spectrum of the 75-foot grating of 
the 150-foot tower telescope on Mount Wilson. 

From accordant measures of the sky on several nights 
the following values were obtained: 


From 1 filter, 10,427A, weight 1; 
From 2 filters, 10,452A, weight 2; 
Mean 10,444A+10A. 


The assigned probable error does not allow for other 
nearby lines or bands, but the main radiation is so strong 
compared with everything else in the whole range from 
3500A to 12,000A that probably no other radiation is 
effective. 

The only possible identifications that suggest themselves 
for this night sky emission are: 


(a) the forbidden *D—*P multiplet of NJ, \10,407.3 
—\10,397.8; 
(b) the (0,0) band of the first positive group of Ne. 


The [NJ] lines do not provide a satisfactory identifica- 
tion. Their mean wave-length 10,403 is definitely shorter 
than that of the observed infra-red radiation. Moreover, 
their total intensity should only be 49 times that of the 
4S—*P line at \3466.4 which is absent from or hardly 
visible in the night sky. Hence the [NJ] lines could not 
account for the intense radiation observed. 

If emitted in the upper atmosphere, the (0,0) band of 
N; should have its intensity maximum somewhere near 
10,450. Failure to observe other strong N: bands indicates 
in the night sky a mechanism enhancing the (0,0) band 
relative to the other vibrational transitions. It is suggested 
that during the night N: molecules are brought in the 
level v’=0 of the B*x, state in three-body recombinations 

The nitrogen atoms would be formed during the day 
through photo-dissociation of N; by solar radiation. The 
suggested mechanism can be considered only for the 
value D(N:)=7.38 v of the heat of dissociation of N: 
which agrees almost exactly with that of the v’=0 level 
of the B*x, electronic state. In fact, the energy D=7.38 v 
is just a little larger than B*x, v'=0, but smaller than 
v' =1. 

This mechanism, which is similar to that suggested by 
Chapman for the [OJ] night sky emission, implies the 
presence of a fairly large number of nitrogen atoms in 
the upper atmosphere. Although the photo-dissociation of 
N: by solar radiation is usually considered as very weak 
compared to that of O; it is by no means negligible. This 
question has been considered recently by Ta-You Wu.' To 
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the arguments in favor of the presence of nitrogen atoms in 
the upper atmosphere advanced by Wu one should add 
that Nz molecules can be predissociated by radiation of 
wave-length shorter than approximately 1400, according 
to a private communication of Dr. Herzberg. 

The mechanism suggested here depends on the value 
D(N:2) =7.38 v; it would not be effective if the heat of 
dissociation D(N:)=9.76 v, recently suggested by A. G. 
Gaydon? were adopted. 

Details will appear in a forthcoming paper in the Astro- 
physical Journal. 

1 
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An Observation on the Ionosphere During the 
Solar Eclipse of July 20, 1944 


Cc. K. Jen, K. T. Cuow, Y. T. Lo, anp C. C. Kuan 


Radio Research Institute, National Tsing Hua University, 
Kunming, China 


September 20, 1944 


N observation of ionospheric changes was made here 
during the solar eclipse of July 20, in the midst of a 
series of routine observations already in progress for some 
time. The eclipse was annular along the belt in southern 
Asia, including such places as India, Burma, and Indo- 
China. As viewed in Kunming, it was only partial, with a 
maximum of 76 percent in diameter. The transmitter and 
the receiver were both placed in a suburb of Kunming, 
separated by a distance of about two kilometers. The 
transmitter is of a conventional build for pulse trans- 
mission, having a peak power of about 500 w and an ad- 
justable frequency range of 2-5 mc. The width of the pulse 
is about 20u and there are 50 such pulses per second. The 
receiver is designed to accommodate a wide-band reception, 
so as to preserve the sharpness of the received pulses. 
The signal and its reflections are observed on the screen of 
an oscilloscope (with adjustable linear sweep and extra 
time calibrating arrangements). For the week preceding 
the eclipse, the transmitter was continuously operating on 


_ the frequency of 3 mc, and for the sake of continuity, it 


was decided to let it stay put on the eclipse day and a 
month thereafter. We have noted that for this frequency 
only the reflections from the E layer have been most 
consistently observed; the reflections from the F layer 
only begin to occur more prominently during late after- 
noons. The virtual height of the EZ layer has always re- 
mained sensibly constant around the value of 120 km, as 
much on the eclipse day as any other day. Since this is 
only a confirmation of the results obtained by many other 
experimenters, we do not think it worth while to show our 
results. In Fig. 1, there is shown the record of the amplitude 
variation of the EZ layer reflection on the eclipse day, as 
compared with that of the average of seven typical records 
taken during the days in the neighborhood of the eclipse. 


THE EDITOR 


Conn 

L 

\ 

o 

7 

=! 
7 10 12 13 14 


KUNMING LOCAL TIME 


Fic. 1. carve for amplitudes of selection on the day 

uly 20); ¥ curve for amplitudes averaged from seven control 

days des of reflections are as fractions relative to 
pow of the ground wave. 


The amplitude of each reflection represents the height to 
which the reflected pulse comes up most frequently (the 
echo wave has a tendency to fluctuate up and down) and 
is oftentimes smaller than a certain maximum peak that 
occurs less frequently. It is seen that on the average the 
amplitude of E reflection generally has a small peak just 
before noon but soon dies down thereafter and only begins 
to rise as the afternoon progresses. The record for the 
eclipse day is considered as remarkable in that it is very 
much different from the usual pattern and its rise and fall 
are fairly coincident with the beginning and end of eclipse, 
The general elevation of the reflected amplitude is most 
extraordinary. There is, further, a dip of amplitude, 
strangely coincident with the maximum eclipse. There was 
absolutely no wishful juggling of the result, because we 
had no idea which way to expect. These results evidently 
suggest a large decrease of absorption with respect to the 
E reflection during the solar eclipse. This diminishing of 
absorption may be caused by the changes occurring in the 
E layer itself as the result of decreased ionization. Or it 
may well be caused by a lower absorption layer that may 
control the amplitude variation. One may suppose that the 
rapid de-ionization and recombination processes that may 
take place in this hypothetical layer may be responsible 
for the decreased absorption. The constancy of the virtual 
height of the EZ layer and the extraordinary increase of the 
reflection amplitude seem to lend help to this hypothesis. 
But the mysterious middle dip is surely not readily explain- 
able on this exceedingly crude model and there must be 
many other difficulties in explaining various other phe- 
nomena. The present experimental evidence is totally in- 
adequate to decide between the two possibilities, but it is 
not impossible that the issue could be settled by studying 
more of the amplitude changes in the future. 
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